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ABSTRACT 


This thesls lays the foundation for the dynamic. stall) 
Investigation belng conducted at the Fluid Mechanics’ Lab- 
oratory at NASA-Ames Research Center. Using existing optical 
and electrical equipment, a new dedicated Mlcro-VAX computer 
and Labstar software, an Indraft transonic wind tunnel = and 
able techniclans to make the proper interface hardware, the 
project came together In a new test facilli!lty at the Fluid 
Mechanics Laboratory. The goal of the thesis was to obtain 
both qualitative and quantitative informatlon about the wake 
profiles of an alrfoll In steady state operations at varyling 
angles of attack and tunnel conditions. To accompl!ish this 
task, sgchlleren photography was used to obtaln a qualitative 
picture of the flow fleld. With thls Informatlon, a two com- 
ponent Laser Doppler Velocimeter was get up to accurately 
measure the velocity proflles that correspond to the schlle- 
ren photographs. Once this preliminary work is completed, 
the same apparatus wll] be used to further Investigate the 


unsteady dynamic stal! phenomenon. 
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I. INTRODUCTION 


A. DESCRIPTION OF THE OVERALL DYNAMIC STALL PROJECT 
Static stal!] Is a well understood and documented phenom- 


enon. A wing or rotor blade that 18S operating at a high a, 


angle of attack, will experlence a point where separated 
flow wll] cause an accelerated decrease in J]Jift anda large 
rapid increase in drag. This ea at which the wing wil) pro- 


duce the maximum /ift, Climax? Is a very Important aspect 
of airfoll performance. This stall] @ determines the stall] 
speed of the alrcraft. Considerable effort has been expended 
to find ways to decrease the stall] speed by increasing the 
Climax: The majority of this effort has been directed to- 
ward a gradual] approach to the stall a, In other words, 
static stall research. 

The wing can be made to continue to produce |!ift past 
the static stall awa lf the wing angle of attack Is In- 
creased rapidly past the stall «a. The phenomenon of 
dynamic stall Is manifested In the rapid decrease of the 
lift at the higher stall a. There are four basic stages of 
dynamic stall] that have been documented (Ref. 1]: 


1. Tralling-edge flow reversal, progressing forward along 


the airfoll. 
2. Formation, growth and shedding of the vortex. 
3. Full] stall. 
4. Boundary layer reattachment. 
The understanding of this phenomenon wil) be beneficial to 


10 


the helicopter, where constantly oscillating rotor blades 
are consistentiy In_ the ae stall region throughout 
their revolution, and Jjet fighters, where agile and rapid 
maneuvering through the stall a t!{s required for surviva- 
bility. By quantifying the velocity field around and behind 
the airfoil, a better understanding of the dynamics of this 
phenomenon will be gained. This new analytical information 
will provide valuable Information to the design engineer to 
expand the flight envelope of future fighter aircraft and 


improve the performance of the helicopter rotor system. 


B. THESIS GOALS 

This thesis lays the foundation for the dynamic” stall 
investigation project being conducted at the Fluid Mechanics 
Laboratory at NASA-Ames Research Center. Using existing 
optical and electrical equipment, a new dedicated Micro-VAX 
computer, new Labstar software, an existing transonic In- 
draft wind tunnel and able technicians to make the proper 
interfaces, the project came together ina new test cel! of 
the Fluld Mechanics Laboratory. The goal of this thesis was 
to first obtain stroboscopic schlieren photographs of a sta- 
tlonary airfoll at different a’s to obtain qualitative 
data on the flow past the airfoll and the size, shape and 
extent of the wake and then set up a Laser Doppler Veloci- 
meter to obtain quantitative velocity data in order to con- 


struct an accurate velocity profile of the flow field. 
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C. SCHLIEREN AND LASER DOPPLER VELOCIMETRY BASICS 
1. Schlieren 

Many optical flow visualization techniques have been 
developed for gas flows. Each of these techniques are sensi- 
tive to variations In density and derlvatives of density In 
the flow. Interferometry !s sensitive to the density of the 
gas flows; ¢. Schileren 1s sensitive to the density grad- 
lent or the change t!n density with respect to distance; 
S$e¢78y. Shadowgraphy Is sensitive to the second derivative 
of the density; S2er $y. All the techniques rely en 
the fact that the speed of Iight !s a functlon of the local 
Index of refraction In that medium. In gases, the Index of 
refraction is a functlon of density. 

The schlieren, German for streak, technique allows 
one to view, qualitatively, the density gradient In a given 
flow field. This visualization system !s considerably sim- 
pler and cheaper to use than other techniques. Therefore, !t 
was chosen as the preliminary technique to obtain qual !ita- 
tive information about the alr flow around and in the wake 
of the alrfoll. The schlleren setup Involves using a Ine 
source of light and projecting It through the test section 
where the flow fleld alters the path of the |ight. The I]Iine 
source |Is usually a monochromatic gas discharge tube. The 
light is projected through a narrow slit to provide a line 
of light of uniform |]lumination. The dimensions of the slit 


control] the definition of the Image seen on the screen. A 


Ww 


high quality concave mirror is usually used to direct the 
light from the source through the — section. On the other 
Side of the test section, another similar mirror is used to 
reflect the light to a viewing screen or camera. Each point 
of the line source will] illuminate each point in the test 
section if the source 1s collimated and of equal optical 
distance on both sides of the test section. The light from 
the second mirror is focused onto a knife edge. The place- 
ment and adjustment of the knife edge !s very critical as 
the sensitivity of the entire system depends on it. The 
knife edge is placed In the plane of the source Image and 
interferes with a portion of this Image. A portion of the 
light from each point In the test section is blocked which 
decreases the light Impinging on the viewing screen. The 
density gradients encountered in the test section normal to 
the knife edge cause a refraction of light allowing more 
light from that portion of the test section to pass over the 
knife edge and illuminate the viewing screen. By rotating 
both the line source and the knife edge both directional 
derivatives of density can be viewed. Also, through the 
lateral adjustment of this knife edge at the source’ image, 
a fine control over the degree of detal!l resulting from this 
technique can be exercised. 

Tne viewing screen can be any image recording 
device. Polaroid cameras can be set up for instantaneous 


results or 35 mm cameras can be mounted at the focal length 
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of the lens to obtaln numerous photographs during the 
testing. Video cameras can also be used. 
2. Laser Doppler Velocimetry 

The techniques for measuring velocity In the past 
used sensors, which need to be Introduced into the_ flow. 
This necessarily meant tIntruding Into and affecting the 
flow. During the last two decades, a new technique, Laser 
Doppler Velocimetry or LDV, has been avallable. LDV’s are 
Indeed a better technique. It Is nonintrusive and is capable 
of directly measuring velocities of particles In the flow. 
Laser Doppler Velocimetry has been very successful In meas- 
ucing flow velocities in turbulent and reverse flows. It has 
been successfully used for the past two decades to obtain 
velocity measurements accurately. 

Laser Doppler Velocimetry uses two coherent! jjight 
beams from the same source, which Intersect to form the 
measurement volume. The two light beams form an interference 
fringe system that Is crossed by particles moving with the 
flow. The light scattered by the particles is Doppler 
shifted“. The frequency of this scattered light Is 


directly proportional to the particle velocity. With 


1. Coherent Light = light In which the electromagnetic 
waves maintain a flxed phase relationship. CRandom House 
College Dictionary) 


2. Doppler Shift = the magnitude of the apparent change 


In frequency resulting from the relative motlon of the 
light source and the recelver. CRandom House College Dict) 
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suitable selection of flow particles, Known also as seeds, 
the fluid velocities can be measured from the particle 
velocities. Using photodetectors, the Iight scattered by the 
particles is converted into electrical signals. This signal, 
after processing, provides frequency information from which 
the velocity can be determined. The measurement of two 
components of the velocity can be obtained by setting up two 
orthogonal systems with laser beams of two colors, In the 
same optics train, thus !dentifying each component with a 


particular color. 
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A. SCHLIEREN SETUP 

The schlieren system for this experiment was set up 
using both a continuous Iight source anda stroboscopic 
source. The continuous source was used to provide a time 
averaged schlieren photograph and also to allow use of a 
video tape recorder to record the time history of the flow. 
To provide instantaneous schlieren photographs of the _ flow 
field, the stroboscope was used. See Figure 2.1 for the 
schematic of the setup. 

To obtain the 35 mm photographs of the flow fleld, a US 
Sclentific Strobrite Stroboscope was used. The strobe had 
five settings, 20 Hz to 1 KHz, to adjust the duration of the 
strobe. The 70 Hz setting proved to be the optimal setting. 
A mechanically adjustable slit in front of the strobe was 
used to produce a |!ine source of |!Iight needed for the tech- 
nique to operate effectively. A large concave mirror with a 
305 cm focal length reflects this lIlight sheet towards the 
test section. The test section had optical access’ through 
both sides of the tunnel. The test sgection’s dimensions were 
13 cm by 25 cm and was constructed out of 2.5 cm thick 
plexiglas. The light passing through the test sgect!ion was 
then reflected by an Identical mirror located opposite the 


first mirror on the other side of the test section as shown 
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Figure 2.1 Schileren Schematic 


imeem 2.2, bolgme from the second mirror was focussed on 
a vertical knife edge whose posit!ion could be adjusted to 
get a highly sensitive system. The knife edge was oriented 
vertically to obtain the visual signature of the density 


gradient in the vertical! plane of the test section. After 
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the knife edge, a 200 mm lens was.used to focus the image 
onto the photographic device. 

Since the test section is 25 cm wide, the photographs 
obtained from the system actually represent density gradient 
effects averaged over the test section width, as is always 
the case with the schlleren technique. 

A Nikon 35 mm camera body with motordrive Cbut no lens) 
was used to record the schlieren Images. The camera was 
placed at the focal polnt of the 200 mm lens. Fine adjust- 
ments to focus the image !n the camera were made by a micro- 
meter slide on the camera mounting mechanism. Pictures were 
taken by opening the shutter with a remote shutter release 
cable, before the strobe was fired remotely and then closing 
the shutter. With a IJlittle practice, this procedure was 
quite rapld. Photographs were taken at two different Mach 
numbers, 0.30 and 0.45, and at three different a’s; 0O9, 
7.79, 15.59. See Figures A.1 through A.6, Appendix A, 
for examples of schlieren photographs. This series of photo- 
graphs was studied to analyze the flow qualitatively and to 
determine the locations In the wake for obtaining the veloc- 


ity data for the two different Mach numbers. 


B. LASER DOPPLER VELOCIMETER SETUP 
This section will discuss the specific LDV system setup 
In test cell] number four at the Fluld Mechanics Laboratory. 


The following sections wlll be discussed: transmitting 
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optics, receiving optics and photodetectors, test section, 
signal processing and data processing. Figure 2.2 shows a 
photograph of the transmitting optics and Figure 2.3 shows a 
photograph of the recelving side. Individual schematics wil! 
be discussed In detall within the following sections. 

The laser used !s a Spectra Physics, Model 164, 5 Watt 
Argon-Ion Laser capable of outputting light from 0.351 to 
0.528 microns and up to twelve discrete wavelengths simul- 
taneously (Ref. 2:p. 23). By choosing the blue beam, 0.4880 
microns, and the green beam, 0.5145 microns, for the two 
component Laser Doppler Velocimetry System, nearly 75% of 
the Initial output power could be used. This’ standard 
practice was adapted in this case. 

1. Transmitting Optics 

Owlng to traversing requirements, the laser was 
higher on the optics table than the 10.8 cm optical! height 
of the transmitting optics. A set of two mirrors at 45° 
angles to the direction of the Incident light were used to 
lower the beam to the require height. Figure 2.4 shows a de- 
talled schematic of the entire transmitting optics system. 

The vertically polar!zed beam exiting the multiline 
laser, after being directed by mirrors, is passed through a 
beam collimator, part of the TSI Mode! 9106 set. Incident 
light Is collimated !f the beam diameter does not change as 
it passes through space. In reality, however, it is suffic- 


lent if the collimator can insure this quality through the 
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Figure 2.2 Transmitting Optics Photograph 





Pnotograpn 


3 Receiving Optics 


Figure 2. 
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Figure 2.4 Transmitting Optics Schematic 


relatively short optical path of the system. Proper colli- 
mation is said to be achieved when the waist! of the beam 
colncides with the focal point of the transmitting lens. By 
Sliding the collimator sleeve in and out and observing the 


beam "spots" on a distant wall, the collimator can be set. 


1. Walst = the polnt on a beam where the diameter is 
the smallest. 


21 


From the coll!imator, the beam passes through a po- 
larization rotator and an attenuator. The polarization ro- 
tator changes the polarization of the Incoming laser beam to 
conform to the polarlzatlon requirements of the next optical 
component, a dispersion prism. The attenuator works In 
tandem with the rotator to insure passage of only the horiz- 
ontally polarized light through these components. The dis- 
persion prism separates the wavelengths, oor colors, of the 
light. The two primary colors of this laser, green (0.5145 
microns) and blue (0.4880 microns), are used to form the two 
component system. A set of mirrors, TSI Model 9107, is used 
to direct each of the beams to the appropriate axis for the 
transmitting optics as shown In Figure 2.4. As can be seen 
from the figure, the green beam is directed down the center- 
line of the optics and the blue !s directed 25 mm to the 
right of the green. 

The centerline, green, beam passes through a polar- 
Ization rotator, TSI Model] 9102, to rotate the beam polar!- 
zation 90°, then a beam splitter, TSI Model 9115. The 
polarization change is necessary to Insure proper splitting 
into two beams from the Incident light. The beam splitter 
divides the Incident beam Into two beams of equal Intensity 


and parallel to the optical axis, spaced 50 mm apart. 
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Since the flow fleld to be Investigated will] have a 
turbulent wake with possible separated flow regions, a freq- 
uency shift system is used for each of the two colored 
beams. Upshifting of the frequency of the laser l]ight prior 
to entering the test section will allow detection of the 
direction the particle traverses through the probe volume. 
Tne shift results In a moving fringe pattern that will gen- 
erate a photodetector signal equal to the shift frequency 
when measuring a particle with zero velocity. This system 
allowS a greater range of measurements of negative veloc- 
ities and, thus, !s required for this system to obtain valid 
measurements. 

The pair of green beams now enter a frequency shift- 
ing device called the Bragg cel!. This device is an acousto- 
optic unlt that upshifts the frequency of an Incident beam. 
The shift is set at 40 MHz. This device is placed i!In_ the 
optical path of one of the two beams forming the probe vol- 
ume for that color. The outgoing beam wil!l1 be shifted by a 
40 MHz shift relative to the other beam. The Bragg cel! 
operates by deflecting one beam in a fan-like pattern In 
the plane intersecting both Incident beams. Each separate 
deflected beam has a frequency that is a multiple of 40 MHz 
from the incident beam. The Bragg cel! should be adjusted 
to shift the light intensity to the first order beam. This 
is accomplished by tilting the Bragg ce)] In tts mount. 


This first order beam should have an intensity greater than 
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90% of the incident light beam (Ref. Sip. 148). An optical 
wedge is included In this device at the output end of the 
Bragg Cell to steer the deflected, first order beam para) le! 
to the unshifted beam. The result of this Bragg shifting Is 
an interference pattern In the probe volume that will move 
relative to a fixed particle. This moving fringe pattern can 
be used to differentiate the direction of motion through the 
probe volume. Particles moving with the direction of motion 
of the fringes gives a lower frequency and against the 
fringe movement gives a higher frequency as shown in _ the 
following equation. The particle moving through the probe 
volume will produce a signal whose frequency lis 


f 


signal = !fgniet +7- fDoppler| 


where the sign !s + for particles moving in the opposite 
direction to the fringe movement (Ref. 2:p. 67]. 

For optic tube support purposes a TSI, Model 9176, 
Support is used. This device also serves to increase the 
effectiveness of thé beam Stops of the blue béam by ie 
reasing the deflected beam separation distance. 

The off-axis blue beam is now brought to the center- 
line using a beam displacer, TSI Model 9174. This model = em- 
ploys a rhombold prism. It displaces the beam a fixed, 25 
mm, distance to the centerline of the optics while keeping 


the beam parallel to the axis. The blue beam then passes 
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through the same series of optical components as the green 
beams: polarization rotator, beam io Litter and Bragg cell. 

The next two elements in the optics train are beam 
steering modules, TSI Mode] 9175, one for each of the 
colors. These modules allow for steering of the shifted beam 
to insure proper beams crossing. 

The Bragg cel!] produces’ several beams, of which 
only the first harmonic, the beam that Is shifted by 40 
MHz, Is important. Al] the other unwanted beams are blocked 
off by means of a beam stop. In the present case, a beam 
stop was used on each of the two shifted beams. 

Following the beam stops, beam spacers, TSI Mode] 
9113-13, are used to for each set of beams to reduce the 
distance between the beam pairs forming the probe’ volume. 
The spacing was reduced to 13mm by using these prisms; 
while maintaining polarity, coherence, beam diameter and 
parallel optical paths. The primary use of these elements is 
to control the fringe spacing In the probe volume by 
changing the angle of intersection of the beam. 
Additionally, however, they permit transmission of the beams 
Into the test section by requiring only sma!! mirrors and 
optic components. 

The four beams from the optic tube are directed to 
the focusing lens with five mirrors. The first two move the 
beams from 10.8 cm to 23 cm above the optical table and the 


next two reflect the beams back 180° toward the last 
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miccor in this series. This mirror directs the four beams 
vertically toward the transmitting optics lens. The four 
beams then enter the transmitting optics lens to focus’ the 
four beams to a point. A final mirror !s used to direct this 
focused beam to the test section. The focusing lens and the 
last three mirrors are mounted ona traversing mechanism to 
enable traversing of the focused beams In any of the three 
orthogonal directions as shown In Figure 2.6. 
<. Receiving Optics and Photodetector 

All the receiving optics are mounted on a traversing 
mechanism which is slaved to the traversing mechanism on the 
transmitting side. The two sides were aligned to within 0.01 
mn and could be moved together. A Taskmaster controller is 
used to move the two traverses both manually and through the 
use of computer programs. See Figure 2.5. 

The receiving optics are used to capture the Iight 
scattered by the particles crossing the probe volume and 
focus this light onto photodetectors. It consists of several! 
components; that are mounted on a hardware arm connected to 
the traversing mechanism. The scattered I!!ight Is collected 
by a series of three lenses that collimate the light and 
direct it toward the polarization color separator. The last 
lens acts as a spatial filter tin addition to being a beam 
collimator. Its function is to allow light focused on the 
lens to pass through while blocking all the rest of the 


light. It Is used to reduce the stray light scattered from 
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around the probe volume, as js the case in measurements 
close to a solld surface. 

The polarization beam separator separates the two 
wavelengths present in the scattered Iight Into its con- 
stituent colors; green and blue. This splitter reflects 
light with horizontal polarization upward while transmitting 
the vertically polarized light through It. This insures sep- 
aration of the two primary colors in the scattered |Jight. 

Two notch fllters are used at this stage, on each of 


the colors to eliminate residual color contamination. The 
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Figure 2.5 Recelving Optics Schematic 
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Separated scattered light 13 now focused on to the photo- 
detector aperture with a lens. 

Tne photodetector converts the Il!ight signal Into = an 
electric signal. The photodetector used !s a photomultipl]ler 
tube or PMT. The advantages of the PMT over the other types 
ace low noise, Nigh gain and broad bandwidth while its 
fragility and the need for a high voltage source are its 
disadvantages (Ref. @2:p. 45]. 

3. Test Section 

The transonic test section is formed by inserting 

upper and lower ventilated walls IJn a Subsonic Indraft wind 


tunnel. The test section’s dimensions are 13 cm by 25 cm. 





Flgure 2.6 Transmitting Slde of Test Section 
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A 7.5 cm chord NACA 0012 alrfoll was mounted between the 
Side walls at the centerline. This airfoll Is attached at 
Quarter chord to a external device to oscillate the airfol| 
through various angles of attack. The side walls of the test 
section are made of 2.5 cm thick plexiglas to allow complete 
and unobstructed optical access to the alrfoll. See Figure 
2.6. The mechanism to mount the model! partially blocks the 
views of the leading edge of the mode!. However, this was 
not a Iilmitation for the present study. 
4. Signal Processing 

The signal from the photodetector must be processed 
to obtain the frequency Information In the scattered I!ight. 
Figure 2.7/7 shows the block diagram of the signal processing 
components. The signal from the PMT !s amplifled by an EIN 
Mode! 403LA, |I!inear amplifier to boost signal strength by 37 
qdB to transmit the signal over a long distance and to be of 
appropriate strength for signal processing. 

A downmixing circult Is used next to decrease the 
effective frequency shift from 40 MHz to around 10 Mhz to 
Improve the counter resolution. The two elements of the 
downmixer are a Tektron!x SG5S03 Levelled Sine Wave Generator 
and a Hewlett-Packard, Model 10534A, Passive Diode Mixer. 
The expected Doppler frequency should be estimated to. prop- 
erly set the mixer frequency. Using the range of Mach num- 
bers expected as -0.2 to +0.5, the free stream velocity 


range can be calculated using the set of equations In 
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Ssectto III A and 209 C as the. ambient total temper- 
ature. The velocity range resulting from this calculation !s 
-~67 to +167 meters per second. With the fringe spacing of 
19.52 microns and the second equation on page 35, the 
estimated Doppler frequencies are -3.4 to +8.6 MHz. Using 
the equation on page 34, the signal and mixing frequency 
should be expected to fall within the range of 36.6 to 48.6 
MHz. By using a starting mixing frequency of 35 MHz, the 


calculated signal frequency reaching the counters will 
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Figure 2.7 Signal and Data Processors 
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be 1.6 to 13.6 MHz. The mixing frequency should be adjusted 
whenever the tunnel conditions change or tin attempts to 
obtaln better counter resolution. 

The next signal processing un!t encountered is the 
series 3000 Macrodyne LDV Processor. This device has several 
functions: variable high and low pass filters, Input and 
output ports, threshold and overload adjustments and _ freq- 
uency counter with fringe passing selection switch. Both the 
low pass filters, to remove the high frequency nolse, and 
the high pass filters, to reduce the low frequency pedestal 
from the signal, are adjustable to provide the optimal fil- 
tering of the signal. Once the signal has been thus condi- 
tioned, a frequency counter {is used to measure the signal 
frequency and convert !t to digital time data that can be 
sent to the data processor. The counter operates on two 
different systems to validate the Doppler signal: a three 
level detection scheme anda S/8 or 10/16 comparison. The 
three level detection scheme shown In Figure 2.8 must be met 
before a signal is validated. A threshold level !s set on 
the counter to be Just above the level of the steady state 
noise. The signal must pass the following gates to be ac- 
cepted: (1) exceeds the threshold setting €2) crosses zero 
and €3) exceeds a negative threshold. To pass this check, 
the selected number of fringe crossings must sequence 
through each check without error. This detection scheme i!s 


notably effective in noise discrimination. Once a signal is 
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recelved and passes the first check, the counter starts to 
measure the Doppler perlod by opening a time gate and stop- 
ping at the switch set number of cycles. This setting of 
cycles Is able to be selected froma switch on the counter 
control panel. Validated data from the switch being In the 
10/16 cycles position tend to be better for strong Doppler 
bursts because the greatest number of cycles measured tends 
to reduce the effects of nolse in the signal. However, this 
setting also restricts the motion of the particle through 


the probe volume and thus’ Introduces artificial biases of 
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the data. The particle must traverse. the probe volume at 
its largest length to Insure paiidated data for 16 fringe 
crossings. For this reason, present experiments were carried 
out with the switch set to the 5/8 position. 
So. Data Processing Hardware 

Once the data Ils received from the Macrodyne signal 
processors, it must be multiplexed and transferred to the 
computer for velocity calculatlons. Tne NASA LDV Computer 
Interface |s now used to control the multiplexing of data, 
sequencing of data, provide a digital display to track 
successful counts of the Doppler frequency and, generally, 
control the flow of data between the signal processor and 
the Micro-VAX computer. Since a two component Laser Doppler 
Velocimeter is used In this experiment, the multiplexer 
sends four channels of data for each valid Doppler burst. 
The first data word contains the time Information, the 
second contains any status information from the signal and 
the last two words contain the frequency data from each 
component. Due to the limited nature of this experlment, 
the first two words of data are not useful In further 
velocity calculations. Therefore, just the frequencies from 
each particle caused Doppler burst are used to calculate the 
respective particle velocities. 

The Micro-VAX computer stores all acquired data ina 
buffer set up by acquisition programs. The next step after 


data acquisitlon Is to calculate the velocitles. These 
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computer programs reduce the massive amounts of raw data 
obtained into meaningful velocity data that can be output in 
tabular form and graphical form. 

The following !s a simplified path followed by the 
Numerous computer programs to reduce the data. At least 1000 
samples from each component are required to obtain a statis- 
tically steady sample frequency. This frequency data must 
then be converted to the true Doppler frequency by the fol- 
lowing equation: 


fqg=fo- fp + fy 


Where fy is the true frequency resulting from the particle 
traversing the probe volume, f, is the signal frequency 
measured by the counter, fp, Is the Bragg shift frequency 
Of 40 MHz and fy is the frequency of the downmixer. The 
downmixer circult permits decreasing of the effective freq- 
uency shift to values within the 1 to 10 MHz range. The true 
frequency data !s then averaged over the samples taken to 
obtain a mean value and a sample standard deviation. Using 
this mean value, for each channel, all samples that are 
within 2 o, two standard deviations, of the mean are saved 
and the rest are discarded. This procedure accounts for a- 
bout 5% reduction of data and will titnsure samples that are 
too low or high will not blas the final results. A new mean 
and standard deviation are calculated for each channel]. Next 


the probe volume dimensions and parameters are calculated 
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with primary emphas!ls on the fringe spacing equation: 

6=A/ (2 * 3in(@/2)) 
where 6 is the fringe spacing in microns, A is the beam 
wavelength in microns and 9 is the angle of intersection 
of the beams. Knowlng the fringe spacing, the velocity can 
be determined using the relation: 

u= fg * 6 

where u is the velocity of the particle in meters per second 
and fq is frequency of the particle measured by the count- 
ers in MHz. 

Once the velocity components of the particle are 
determined streamwise and in the transverse directions, the 
computer stores the results in a data table for further 
analysis. The probe position is also saved in this data 
fille. The entire process is now repeated at different 


locations of interest In the flow. 


C. LABORATORY APPARATUS 

Besides the Laser Doppler velocimeter and schlieren 
systems described earlier, additional Instrumentation was 
used to obtain the wind tunnel pressure and temperature 
data. A 24 port scanivalve pressure transducer system was 
set up to measure the total and static pressures at various 
points Inside the operating tunnel. This scanivalve system 


was driven by a digital to analog interface, built on the 
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premises, which allowed computer control of the “stepp!ng" 
and "homing" of the scanlivalve. 

Using the Digital to Analog Interface on the MicroVAxX 
computer, these signals to step the scanivalve system were 
transferred from a computer program to the physical unlit for 
action. Once the stepper motor was sent the appropriate slgq- 
nal to move, the Analog to Digital Interface board on the 
computer sampled the pressure transducer signals and loaded 
the data Into the computer for conversion Into pressure data 
and later analysis. 

The atmospheric pressure was recorded from a mercury 
barometer with micrometer scale mounted on the laboratory 
wall near the wind tunnel. The total temperature was record- 
ed on a mercury thermometer mounted on the barometer. Both 
of these measurements wlll soon be replaced with Paro- 
sclentiflc pressure sensors and thermocouples that wil! 
nave a direct Iink wlth the computer. 

A TSI Six-Jet Atomizer, Model 9306, was used to provide 
particle seeding for the LDV system. 

Figure 2.9 shows the general J|aboratory setup and 


location of the data processing and anc!!lary equipment. 


D. DATA PROCESSING 
1. Introduction 
The Laser Doppler velocimeter yields a large amount 


of raw data during Its operation. This data must be stored 
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Figure 2.9 Laboratory Setup 
in the computer in a logical sequence in order to properly 
calculate the velocities and al} other associated data. 
2. Computer Proarams 
The following discussion pertaining to the func- 
tioning of each of the computer programs that was written to 
accomp!ish the task of managing the vast amount of data In 


the most efficient and timely manner. 
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a. Mastee Program 

MASTER controls the overall] flow of data during 
each of the subroutine calls, acquires al) the data and 
reduces it to a usable form. 

First, the wind tunnel conditions are sgampled by 
the use of various subroutines: DA11JiLAB, SETUP, GETDATA, 
DETACH and OUTPUTD. The traversing mechanism is then posi- 
tioned with subroutine HWTRVT to move the probe volume of 
the LDV to the proper location in the flow. Next, the LDV 
data Is acquired with subroutine HWWA. The samples are 
taken, velocities calculated using subroutine VELOCITY and a 
statistical mean and standard deviation are calculated. The 
data Is then reduced with subroutine REDUCE_DATA to preclude 
the use of any velocities greater or less than two standard 
deviations away from the mean. Each of the two velocity com- 
ponents are then recorded along with the position ine the 
flow. The velocity proflles are obtalned from the data using 
post-acquisition processing programs. An _ option is given 
during the acquisition process to view a histogram of the 
resulting velocities to check on the distribution of data 
samples with subroutine HISTOGRAM. A sample histogram from a 
test run is shown in Figure 2.10. 

b. Scanlivalve Programs 

First a signal is’ gent to the scanivalve to 

"home" to port O for the first pressure reading. This 1s 


accomp!!shed by using the subroutine DAI1!J1LAB. Next, a cal] 


38 


VERTICAL VELOCITY HISTOGRAY 


BO 200 
To. OO 
70200 
6o-u0 
60.900 
33. OO 
So. 00 
45.00 
40.00 
55:00 
30°. 09 
251, 00 
ZOnCe 
15 OO 
PO Ue 


PERCENTAGE OF PTS 





eon omer so. fe 6 -6.493ee. 86 12.15 24.44 30.73 40.02 


VELOCITT Ices (mos? 


Figure 2.10 Sample Histogram 


to subroutine SETUP is made to connect the computer’s in- 
ternal timing clock and Analog to Digital Interface board to 
program contro!. Also, the sample acquisition rate and num- 
ber of samples are interactively entered. 

Now, the computer is ready to acquire data from 
the scanivalve system. Since there are 24 ports to sample on 
the scanivalve, a loop Is entered in which samples are 
taken, scanivalve stepped and data reduced during each |oop. 
soubroutines GETDATA and DA11Ji1LAB acquire the data and step 
the scanivalve respectively. DA11J1LAB has an adjustable 
time delay to give time for the pressure in each port to 


equalize after scanivalve stepping. 
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After all 24 ports are sampled, the mean and 
Standard deviation voltages are calculated for each port. 
The calibration pressure attached to port Q and the atmos- 
pheric pressure at port 1 are now entered interactively In 
consistent units; inches of mercury in this case. With these 
two pressures entered, a calibration |line Is calculated and 
the remaining voltages from the scanivalve system are con- 
verted to pressures. Using port 2, settling chamber stagna- 
tion pressure, and port 7, test section static pressure, a 
tunnel Mach number is calculated. See section III A for the 
equations used in these calculations. Next, the total 
temperature is entered and the free stream velocity is 
calculated. 

For use in labelling a program generated plot- 
ting fille for post-acquisition analysis, the angle of attack 
of the test section airfoil and current time are entered. 

The last part of the scanivalve data acquisition 
programs close the loop on the analog to digital interface 
by detaching the clock and driver. A diagnostic output file 
is also generated for post-acquisition analysis. 

c. LDV Data Acquisition Program 

Both the transmitting and receiving optics of 
the Laser Doppler Velocimeter are mounted on three axis 
motorized traverses that are slaved together to move the 
probe volume to any point of interest In the flow. Using a 


pre-existing software driver for the Taskmaster traverse 
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control unit, the optics could be moved at will from the 
computer workstation. 

Although a pre-existing software driver was 
available for the NASA-LDV Computer Interface, significant 
changes had to be made to the code for it to work on the 
present system. The malin Laser Doppler Velocimetry 
acquisition subroutine "“DAMUX" was responsible for loading 
the incoming data from the multiplexer into data arrays. 
Once the raw frequency data is In arrays, it has to _ be 
decoded to correspond to usable numbers. 

The time word output, as 14 bits, by the counter 
has a mantissa and exponent format: with a 10 bit mantissa 
and a 4 bit exponent. The computer reads it as a 16 bit word 
and adds two bits that are not used for data handling. These 
two bits are only used during direct computer links. This 16 
bit word is output as a complementary word that Is reversed; 
that is, the zeros are output as ones and the ones as zeros 
then the entire word is reversed (i.e. the least significant 
bit which is on the right most side bit is moved to the left 
most bit In the word). First, the word is made into its com- 
plement then both the mantissa and the exponent are reversed 
prior to decoding. The exponent is preset from the counter 
front panel switch. The following is the standard 16 bit 


word composition and what actually results from the output: 
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Correct output. 
OLTTeTCuonl. Orr er = 465x27 = 59,520 


Mantissa Exp NU 


Actual output 
O1 o001 0111010001 


NU Exp Mantissa 


Translated by the computer 


0100010111 0100 O01 = 279x244 = 4,464 


First step: Take its complement 
10 1110 1000101110 


NU Exp Mantissa 


Next: reverse the order 

0111010001 O111 O1 = 59,520 

Mantissa Exp NU 
where NU stands for Not Used. In order to fill the sixteen 
bit word, the 14th and 15th bit are added to control data 
flow but have no significance to the value of the data. From 
this last 16 bit word Is calculated the Doppler shift freq- 
uency. The output from the NASA-LDV Computer Interface I!s 
actually representative of the time elapsed for the particle 
to traverse 8 fringes of the probe volume. The Doppler 


perlod Is expressed by the following equation: 
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ty = 1/32 (DoDg...Dp> x 2EXP ns 


where D9..Dy represents the binary mantissa and EXP is 
the preselected exponent set on the counter. From the 
period, the Doppler frequency is calculated using the 
following equation: 


fq = 1000/7t, MHz 


This Doppler frequency is then be used to calculate the 
velocity of the traversing particle. Since frequency 
shifting is used to discriminate the actual direction of the 
particle through the probe volume, the Doppler frequency 
must be adjusted to account for the Bragg and mixing 
frequencies. See equation on page 34 for this relationship. 
Once the actual frequency of the particle is determined, the 
fringe spacing !s calculated by the first equation on page 
35 and the velocity is found using the second equation on 


the same page. 
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d. LDV Support Programs. 

Three subroutines are used to provide the 
supporting functions to the LDV data acquisition. These are 
VELOCITY, HISTOGRAM and REDUCE_DATA. 

VELOCITY is responsible for calculating the 
fringe spacing, extracting the particle frequency from the 
signal and then calculating the two velocities. 

HISTOGRAM takes the velocity data and plots 
a histogram of the percentage of velocities that fall within 
a given range. The x axis of the histogram is divided 
equally into 8 segments from the minimum to maximum velocity 
values. The option to plot a velocity histogram is given 
after each LDV acquisition to spot check the distribution of 
the data. A nearly Gaussian or normal distribution is 
expected if the system is operating properly. This will also 
show the effects of nolse as evidenced by the broad band of 
the histogram values. 

REDUCE_DATA reduces the velocity data to 
only those values falling within two standard deviations of 
the mean. This results ina 4.5% rejection of data which 
will reduce the effect of high frequency nolse on the 
counters and tend to reduce the effect of very large 


particles moving slowly through the probe volume. 
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Ill. DATA ACQUISITION 


A. TUNNEL CONDITIONS 

In order to obtain good correlation between the schlie- 
ren photographs and the Laser Doppler Velocimeter obtained 
data, the operating tunnel conditions would have to _ be 
nearly identical. Two sets of conditions were chosen: one 
corresponding to an approximate Mach number of 0.3 and the 
other 0.5. Once the tunnel was running at operating 
conditions, the computer program MASTER calculates the test 
section Mach number and free stream velocity according to 
the following relationships. Assuming isentropic flow in the 
test section, the section Mach number, M, Is calculated from 
the stagnation pressure, P,, and the static pressure, P, 
using the following relationship: 


M = (CCPL/PIS I-17 Y)-1)xcarty-1)))1%e 


where y = 1.4. The resulting calculated Mach numbers for 
the two conditions are as follows: 0.295 and 0.450. From 
these Mach numbers, the corresponding free stream velocities 


are calculated. 


where a is the sonic velocity. 


a = CyRT)17%2 
where R Is the universal gas constant and T is the static 
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temperature In the test section. The static temperature can 
be calculated from the total temperature, To, measured on 
the wall near the tunnel by the next relationship: 


T= TS/(1 + (Cy-1)/2) *M2) 


B. SCHLIEREN PHOTOGRAPHS 
1. Introduction 

Once the wind tunnel] was running at the first of 
two operating conditions, manometer board readings were 
recorded and Mach number calculated using the _ first 
equation on the preljous page. The schlieren photographs were 
then taken at four different angles of attack and three 
different strobe Intensities to obtaln the optimal 
photographs. 

2. Test Conditions and Procedure 

Using the procedures detailed in section II A, the 
schlieren photographs were taken at the settings and 
conditions |]!sted in Table I, below. 

Upon completion of the wind tunnel runs, the 35mm 
film was developed and printed. From these schlieren 
photographs, the conditions resulting In the best depiction 
of the wake profiles were chosen to represent each Mach 
number and angle of attack. These are annotated with an 
asterisk, *, in Table I, below. 

Using these photographs, a short HP-41 CX calculator 


program was written to relate photographic dimensions to 
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TABLE I. SCHLIEREN PHOTOGRAPH PARAMETERS 


Strobe 
M a (deg) Intensity c¢Hz) 
0.30 0.0 20 
0.30 0.0 70 
0230 0.0 150 
0.30 Va 20 
0.30 Facial 70 
es 0.30 Gal d 150 
* 0.30 152.5 20 
0.38 rSo.o 70 
0.30 £55 150 
0.30 Zocc 20 
0.30 23.6 70 
0.30 23.6 150 
» 0.45 0.0 20 
0.45 OF; 0 70 
0.45 0.0 150 
« 0.45 1 20 
0.45 Td 70 
0.45 im. i150 
0.45 fSa5 20 
& 0.45 (fSseS 70 
0.45 © 150 


real test section dimensions. Once the real dimensions of 
the wake were Known, the Taskmaster traverse mechanism was 
calibrated. This yielded a 1:1 relationship in the chordwise 
direction and a 1:4 relationship in the vertical direction. 
This means that the Taskmaster digital control panel reads 
true values for the X direction, chordwise, and four times 
the real value in the 2 direction, vertical. Arbitrarily 
choosing the tralling edge of the airfoll at 0° a as an 
Initiallzation point for the Taskmaster, a table of LDV 
sampling points was made for the first three angles of 
attack. See Table II, below, for a sample of the Taskmaster 


table at a = 7.79 
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Lables ii. 


TASKMASTER LDV SAMPLING LOCATIONS 





Alpha = 7.7 deg 
2 2.00 7? 6.5°6°5:5 5 wt. 5 de S75 sees le =.) ~1 491.9 = Caen 
205.50 7? 6.50865) 5.0. 58 Ach Oe geome eS -.2 -1 -195 -2"tecao 
20 1.00 7 6.55GH5.5 S YS dies Gecmeero leo =.0 ~1 71.5 -2 teeee 
OO 0 ¢ 6.556655 3 8.5 $3 Seo ees io -.2 ~! -1.5 -20eee 
20 ...00 ? 6.5 Ge5S.5 594.5. 4+ 3155 cme lad -.9 -1 -l.5 920 
19 =-.00 ¢ 6.S9i6eS.0 9 4. lsmeeco mec) Ps ae “1.5 -2 same 
17 ={.00 7 G6. SiGgposon Sa 4 Suede os cmmcmEe. © leo -1.5 -2 72am 
17 -1.50 7? 6.5 605.5 5 4eoe*h 305 7-5 lS -1.5 -2 =a 
11 -2.2 7 6.5 i6685.5 S 4eSes4e SiG cme “200 
25-3200 7 6.5) 6, sygtO SW cuits 5 Bomeceo =é.0 


176 


(Total locations sampled) 


C. LASER DOPPLER VELOCIMETRY DATA ACQUISITION 


1. Calibration of the LDV System 


the Doppler frequency and 


The relationship between 
the velocity of the particle traversing the probe volume 1s 


given by the second equation on page 35. The particle 


velocity is directly proportional to the Doppler frequency 


by a factor of the fringe spacing. By the first equation on 


page 35 , fringe spacing !s only ae function of the 


transmitting lens angle and wavelength of the light. Since 


the wavelength of the light is accurately known to within 
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0.05 microns, or 0.01%, an accurate lens angle !s all] that 
is needed for calibration (Ref. 4]. 

The following relationship was used to calculate 
the lens angle: 


@=2x tan! (Ee dy (2 * £)) 


where E is the beam expansion ratio, 1.0 in this case: d is 
the incident beam spacing, preset at 13 mm with TSI Mode] 
9113-13; and f is the wavelength of the light. 

The lens angle in this case was not accurate 
because the focal length of the focusing lens was not 
known. Therefore, calibratlon was necessary to Insure 
accurate velocity measurements were obtalned. 

The calibration method chosen was to operate the 
Wind tunnel at known conditions, calculate the free stream 
velocity from the scanivalve pressure data and then move 
the probe volume via the Taskmaster to colncide with a 
position in which the free stream velocity could be 
measured. From this procedure, a constant could be applied 
to the velocity equation as follows: 


u = fy *% 6 *%# CCalibration Constant) 


Tnis calibration constant will compensate for the estimated 


focusing lens focal length. 


2. Test Conditions and Procedure 
Once calibration Is completed, the tunnel Is 


brought to the test conditions, the Taskmaster is moved to 
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the first sample polnt and the MASTER computer program is 
run with the adjusted fringe spacing to reflect the 
calibration. 

At each sampling point, 1000 samples will be taken 
to insure statistical accuracy of the mean flow calcula- 
tions. The scheduled sampling points can be geen In Table 
II, page 48. After each Taskmaster traverse, the oscillo- 
scope is checked to Insure no drastic loss of the signal 
occurred during the movement. If the Doppler signal de- 
creases, slight adjJustments to the receiving optics traverse 
may be needed to regain the signa! strength. To optimize the 
data received, slight adJustments may have to be made, also, 
to the selected mixing frequency, mixing frequency amp)i- 
tude, counter threshold setting or high and low pass filter 
Sex tings. 

After Table II data points are sampled, the airfoil! 
angle of attack or test Mach number ‘is changed and the above 


procedure is repeated. 
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IV. RESULTS AND DISCUSSION 


A. SCHLIEREN 

The first thesis goal was attained early enough to 
provide qualitative Information concerning the airflow 
around and in the wake of the alrfoll. Schlieren 35mm 
photographs were produced for both selected mach numbers, 
0.3 and 0.45, and for four different angles of attack; 0°, 
7.72, 15.59 and 23.69%. The best photographs for. each 
Mach number and the first three angles of attack are dis- 
played !n Appendix A, pages 55-57. From these photographs, 
LDV sampling points were chosen as seen from Table II on 
page 48. 

Three turbulent areas are evident on all the schlieren 
photographs. As seen in Figure A.4 on page 56, the airflow 
Is slightly disturbed upstream of the alrfoll due to the 
presence of a pitot-static tube extending Into the alrflow 
Just forward of the airfoll. This tube has been subsequently 
removed to clean up the flow. At the top of the same _ photo- 
graph, a wedge shaped disturbance can be seen. Th!s and the 
lower floor turbulent layer are results of the transonic 
wall slottings found on both the roof and the floor of the 
test section. Although a concern for later unsteady flow 
Studies, these two areas do not Interfere with the wake 


profiles resulting from a steady flow analysis. 
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B. LASER DOPPLER VELOCIMETER 

Due to equipment problems and time constraints, the LDV 
data acquisition part of this thesis was not able to be 
accomplished. All systems are now operational, to include 
computer programs, however, there is inadequate time to 
acquire data and analyze the velocity profiles obtained. 

Calibration runs were made with a limited degree of 
success when the wind tunnel was run under Mach 0.3 and the 
Laser was not powered over 2 watts. The problem with the 
Laser is easily repaired, however, time consuming. The 
upper limit on the wind tunnel Mach number is caused by the 
evaporation of the seeding particles in high speed flow. 
This, too, can be easily overcome by the proper selection 


of seeding particles and related atomizer. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The thesis provided invaluable assistance to the dynamic 
Stal] project at the Fluid Mechanics Laboratory. The schlie- 
cen photography was able to provide the qualitative data on 
the steady state flow field. With this accomp!ished, the 
schlieren system in this |]aboratory was able to be easily 
tailored to the next step: unsteady flow analysis. 

Tne LDV work provided the laboratory with the essential 
computer programs, hardware interfaces and drivers, optics 
systems and auxiliary required equipment to enable the pro- 
ject to proceed with very little additional work to achieve 
the steady state velocity profiles. The ground work has been 
accomplished to enable the project to continue forward with 
the systems on hand to study the effects of dynamic stal]] in 


both the steady and unsteady environments. 


B. RECOMMENDATIONS 

The schlieren system is fully functional, however, the 
Laser Doppler Velocimeter system stil! requires a little 
additional effort to make it productive. The Laser doughnut 


condition needs to be addressed to enable higher powers to 
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be utilized. The types of seeding available to the system 
needs to be addressed to insure operability of the tunne} 
and data acquisition systems above Mach 0.3. With these 
minor fixes anda little fine adjustments of the signal! 


processing equipment, the LDV is ready for action. 
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APPENDIX A: SCHLIEREN PHOTOGRAPHS 


Figure A.1 Mach Number 0.3 
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Figure A.2 Mach Number = 0.45, a = 0° 
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Figure A.3 Mach Number = 0.3, as 7,79 





Figure A.4 Mach Number = 0.45, a = 7.79 
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Figure A.6 Mach Number = 0.45, a@ = 15.59 
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APPENDIX B: COMPUTER PROGRAMS 


PROGRAM MASTER 





Master Progrom written by Maj Richord R Ryles, NPS Student 28 Apr 87 


Last Update 6 Sep 87 


This progrom will gather dota from the scanivalve system, move the 
Taskmaster and ocquire LDV velocity data. 

This progrom should be used with subroutine SETUP, GETDATA, DETACH, OUTPUTD 
HWTRVT, HWWA, HISTOGRAM and VELOCITY. 





INCLUDE 'SYS$LIBRARY:LIOSET.FOR’ 





Declare variables 





INTEGER DEV_IO,CLK_ID 
INTEGER NCHANL,HISTO,8LUE_GREEN 

INTEGER NPT,NPTS,NSAMPLES , COUNT 

INTEGER 1,J,NSTEPS, IER, ANSWER, ENO, TALLY, TIME 

REAL VOLTAGES (500@) ,U_U( 260) , V_V(20@) ,U_V(20@) ,DELX ,DELZ 

REAL®8 SUM,AMEAN,X, X,Y, YY,STD_DEV,.STO_DEV_V,STO_DEV_U, SUM_V, SUM_U 
REAL*8 SUM_VERT,SUM_HORIZ ~ 

REAL RATE, SCANI_VOLT(@: 23,2), SCANI_PRESS (0:23) , VERT_VEL (8200) 
REAL CALIB_PRESS, SLOPE, INTERCEPT ,HORIZ_VEL(8200) , FINAL_V,FINAL_U 
REAL X_POSITION, Z_POSI TION, MEAN(2) ,.M, FREE_STREAM, TO, ATMOS_PRESS 
REAL X_LOC(20@) , Z_LOC (208) ,U( 206) .V(20@) , ALPHA, ZERO, CORR_FACTOR 
REAL*&® UU,VW,UV 

CHARACTERe1 ANS, VOLT_ANS 

LOGICAL HOME 





Scanivalve data acquisition | 





OPEN (UNI T#99,FILE=’ RESULTS .DAT’ ,STATUS=’ NEW’ ) 
TALLY = 0 

HOME = . TRUE. 

NSTEPS = @ 

CALL DA1IJ1LAB (HOME ,NSTEPS, IER) 

IF (IER.NE.1) GOTO 178 

ANSWER = 6 

WRITE (¢,°) ‘Skip to Taskmaster? (y/n) [DEF - No]’ 
READ (¢,1) ANS 

IF (ANS.EQ.'Y’.OR.ANS.EQ.’y’) GOTO 90 

WRITE (e,¢) ‘Print the scanivalve diagnostic data to a file? (Y/N)’, 


& ' [DEF - Mo)’ 


READ («,1) VOLT_ANS 
CALL SETUP (CLK_IO,OEV_IO,RATE,NCHANL NPT .NPTS) 
CONT INUE 
HOME = .FALSE. 
NSTEPS = 1 
DO 50 J = 6, 23 
WRITE es J 
IF aera CALL DATIJILAB (HOME ,.NSTEPS, IER) 
IF (IER.NE.1) GOTO 158 
CONT INUE 
CALL GETOATA (OEV_ID,NPT ,NPTS, VOLTAGES) 
*®* @eeeeeeeeeeeeee#ee#éeée@ 
Calculate mean and stondord deviation 


e 


Ceeeveev4nsresvseeeeveeeeneseee @ 


SUM = 6.0 
00 3@ 1 = 2, NPT 
SUM = SUM + VOLTAGES(I) 


08 
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aQa0 
. 
e 


Es(1) — AMEAN 
e (1.0/FLOAT(NPT-1)) 


~< 


STO_DEV = SORT (SUM) 
IF (VOLT_ANS.EQ.'y’.OR.VOLT_ANS.EQ.‘Y’) CALL OUTPUTO (RATE,NPT, 
NCHANL , AMEAN, STO_DEV, J, VOLTAGES ) 

Sa an es = AMEAN 

SCANI_VOLT(J,2) = STO_DEV 
CALL DETACH (CLK_IO,OEV_ID) 
IF(VOLT_ANS.EQ.’Y’.OR.VOLT_ANS.EQ.’y') THEN 

OPEN(UNIT=15, FILE=’SCANI_VOLT .DAT’ ,STATUS#='NEW’ ) 

DO 66 I = 8, 23 

WRITE (15,3) I, SCANI_VOLT(I,1),. SCANI_VOLT(I,2) 
ENO IF 
eeee#sese#e#e#eeeeee%ee#ee#eee#e#eseeee#es@é@ 
Calculate scanivalve calibration line 
and convert all voltages to pressure 
e®eeee#e#28ke#eee8e86dc aMmUmUCUHMDmUCUCOUCOHDUC SUC‘ SUC SIU UIUCUCOULUCUDClCU 
WRITE (¢,¢) ‘Enter Scanivalve Calibration Pressure: ' 
READ ee CALIB_PRESS 
WRITE(e,*)’Enter Calibration Pressure Correction Faoctor:' 
READ(«,*)CORR_FACTOR 
CALIB_PRESS = CALIB_PRESS + CORR_FACTOR 
WRITE(e,*)’Enter atmospheric pressure: 
READ(e,¢) ATMOS_PRESS 
SLOPE = (ATMOS_PRESS—CALIB_PRESS )/(SCANI_VOLT(1,1)-SCANI_VOLT(@,1)) 
INTERCEPT = CALIB_PRESS - SCANI_VOLT(®,1)*SLOPE 
WRITE(99,¢)’CALIB_PRESS’ ,CALIB_PRESS 
WRITE(99,¢)’CORR_FACTOR’ ,CORR_FACTOR 
WRITE(99,¢)’ATMOS_PRESS’ ,ATMOS_PRESS 
OPEN(UNIT=16,FILE=’SCANIPORT.DAT’ ,STATUS='NEW' ) 
DO 70 I = 6, 23 
SCANI_PRESS (1) = SCANI_VOLT(I,1) © SLOPE + INTERCEPT 

WRITE(16,¢)1, SCANI_PRESS(1) 
CLOSE( 16) 
IF (VOLT_ANS.EQ.’°Y’ .OR.VOLT_ANS.EQ.'y’) THEN 

DO 88 I = @, 23 

WRITE (15,4) I. SCANI_PRESS(1) ,SCANI_VOLT(I,2)eSLOPE+INTERCEPT 


END IF 

CLOSE(15) 

WRITE(e,¢)° ° 

WRITE £2) “equat ion of the calibration line =’ 
WRITE(e,¢)'Y @’ slope,’ X +’, intercept 


WRITE(99,¢)'Y @’, slope,’ X +°,intercept 

WRITE(*,*)° ° 

eoe#eeee#2#ee@e@mm.em6mUmUOmWmUUC OWUC HMUCDUC“ SUC SLU UU OUCOUCUCOUCUCOUClCU 

Calculate tunnel Moch and velocity 
oeoees8e?8Fe2e#e8080h6hm8mCmtUCUHMUC HUCHUC HDLC OHUCMOOrmhUCOrhUC HTWhCU(CU ]mhUlUhHhlUw! 

IF (SCANI_PRESS(7).E£Q.@.0) SCANI_PRESS(7) = 1.2 

IF (SCANI_PRESS(2).LT.SCANI_PRESS(7)) GOTO 85 

M = SORT(((SCANI_PRESS (2) /SCANI_PRESS(7) )ee(2./7. )-1.)#5.) 
WRITE(e¢,*)’Tunnel Moch Number @’ ,M 

WRITE(99,¢)'Tunnel Moch Number =’ ,M 


are neni ; 

WRITE(e,*)’Enter Total Temperature (C)' 
READ(e,¢) TO 

WRITE(99,¢)’° TOTAL TEMP’, TO 
WRITE(*#,0)* ° 


TO = TO + 273.15 

FREE_STREAM = M ¢ SORT(401.520T0/(1.+(Mee2/5. ))) 

ee PE tsoRanurte Streom Velocity (mps) =’, FREE_STREAM 
WRITE(99,¢)’Free Streom Velocity (mps) =’, FREE_STREAM 
IF(FREE_STREAM.EQ.0.0) FREE_STREAM = 1.£-10 
WRITE(e,¢)’ ° 
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OOOO 


100 


aa 
“ 


e@d ee e 


NAAQAAAQA=0O0OAN 
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WRITE(e,¢)’Enter angle of attack (deg): ’ 
READ(«,¢) ALPHA 

mil telesty Gat eme ALPHA 

WRITE(e,*)’Enter Current Time: (4 digit integer)’ 
Read(e,¢) TIME 

SI TER oe cee 

WRITE (e,¢) ‘Skip to end? (y/n) [DEF - No)’ 

READ (¢,1) ANS 

IF (ANS.EQ.'y’ .OR.ANS.EQ.'Y’) GOTO 168 


Taskmoster mover and LOV data acquisition | 





WRITE (¢,¢) ‘Which color is connected to the left ($3) channei?’, 
& 


’ (Blue=! Greens)’ 

READ eoiee BLUE_GREEN 
WRITE(99,¢) 'BLUE_GREEN’ ,BLUE_GREEN 
IF(FREE_STREAM.EQ.0.0) FREE_STREAM = 1.€-6 
OPEN (UNIT#=25,FILE=’FINAL.DAT’ ,STATUS='NEW’ ) 
WRITE Be ? Taskmaster (Final .Dat)’ 
WRITE (25,¢) ’ pos (in) velocity/std dev velocity/std dev’ 
WRITE (25,16) 
WRITE (e,¢) 'Skip to LOV? (y/n) [Def - NO]}’ 
READ (e,1) ANS 
IF (ANS.EQ.’Y’ .OR.ANS.EQ.*y') GOTO 110 
IF(TALLY.EQ.@) THEN 

WRITE(e,¢)’ Enter X—-Position for calibration’ 

READ(e,¢) X_POSITION 

WRITE(e,¢)’ Enter Z-Position for calibration’ 

READ(e,*) Z_POSITION - 
ENOI F 
CALL HWTRVT (DELX,DELZ) 
WRITE (e,¢) ‘Move the Taskmaster again prior to LDV Acquisition? (Y/N) 
READ (¢,1) ANS 
eeeee#e##ee$g#%eeeeee#eeeee#esesestst e@ 
Record position of LOV at all times 
eeeeeeee%ee%ee#e#esesee#eseeeesee?stst ese 
IF(ANS.EQ.’°Y°.OR.ANS.EQ.’y') THEN 

X_POSITION = X_POSITION + DELX 

Z_POSITION = Z_POSITION + DELZ 
GOTO 108 

ENDIF 
@eeeee*eees#e#e¢etc_oeee#ee#eee%#ee%eeeeesste e® 
Acquire LOV data 
eoeeee%e#8e#e##ee%#ee#e#esetsesses#stesteeeeesd @ 
CALL HWWA (NSAMPLES , ANSWER , BLUE_GREEN) 
eoeeesee#eeeeeee#ee#e#e#e#e%seestesresteteeses se 
Increment position and record it 
Read velocity dota, sum it and find 

standard deviation 
e®eeeee#eee#eee@¢egee#eeee#e#ee%eeeee#sest ® 
X_POSITION = X_POSITION + DELX 
Z_POSI TION = Z_POSITION + DELZ 
SUM_VERT = 0.00 
SUM_HORIZ = 6.00 
OPEN (UNIT#35,FILE='’ VELOCITY .DAT’ ,STATUS="OLD’ ) 
00 120 I = 1, NSAMPLES 

IF be aE omen ea am ren ts, NCATE AML CaEa eet 

IF (BLUE_GREEN.EQ.1) READ(35,7) HORIZ_VEL( 1), VERT_VEL(I 

SUM_VERT = SUM_VERT + VERT_VEL(1I) 

SUM_HORIZ = SUM_HORIZ + HORIZ_VEL(I1) 
CLOSE( 35) 
FINAL_V «= SUM VERT / FLOAT (NSAMPLES) 
FINAL_U » SUM_HORIZ / FLOAT(NSAMPLES) 


.0® 
.0® 
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1 


aq 


aan 


3® 


210 


00 13@ I @ 1, NSAMPLES 
X = VERT_VEL(I) - FINALV 
XX = HORIZ_VEL(I) - FINAL_U 
Y =X eX @ (1.0/FLOAT (NSAMPLES—1) ) 
YY = XX © XX © (1.6/FLOAT (NSAMPLES-1)) 
UJ = W + ((HORIZ_VEL(I )ee2)—(FINAL_Uee2)) 
VW = W + | VERT_VEL(1)ee2)-(FINAL_Vee2) 
UV = UV + ((HORIZ_VEL(1)*VERT_VEL(1))—(FINAL_UeFINAL_Y) ) 
SUM_V = Y + SUN_V 
SUM_U = YY + SUM_U 

UU = WU / FLOAT(NSAMPLES 

WV = W / FLOAT(NSAMPLES 

UV = UV / FLOAT(NSAMPLES 


STD_DEV_V = SQRT aU; 
STO_DEV_U = SORT (SUM_U 
WRITE(e,*) ° ' 
WRITE (¢,¢)'Velocities and std devs =’ 
WRITE s3}'u dir’, FINAL_U,STO_DEV_U 
WRITE (e,¢)’V dir’, FINAL_V, STO_DEV_V 
WRITE (99,¢)'Vetocities and std devs =’ 
WRITE 99,0) dir’ ,FINAL_U, STO_DEV_U 
WRITE (99,¢)’V dir’ ,FINAL_V,STO_DEV_V 
WRITE (e,¢)’ ' 
WRITE(e,¢)’First cut RMS values’ 
WRITE 3) "Uo Yor WW 
@,¢)°UV—-RMS’ ,UV 
WRITE(99,¢)’First cut RMS values’ 

3) Go Ru WW 

, ©) UV-RMS’ , UV 
e®*eeeeeeereeee%ereeeeeeeeeeeeeeees ee 
Eliminete doto outside two sigma from mean 
ee oeoeeeeee#eeeeee%eseeseeesveeeees?®@ 
IF (BLUE_GREEN.EQ.1) THEN 

eas = FINAL_U 


MEAN(2) = FINAL_V 
CALL REDUCE_DATA (NSAMPLES , COUNT, ENO ,MEAN,STO_DEV_U, STO_DEV_V, 
BLUE_GREEN) 
FINAL_U = ait 
FINAL_V = MEAN(2 
WRITE (99,0) U air velocities ond std devs =’ 
WRITE (99,¢)°U dir’ ,FINAL_U,STD_DEV_U 
WRITE (99,¢)°V dir’, FINAL_V, STO_DEV_V 
ELSE 
mets = FINAL_U 
MEAN(1) = FINAL_V 
CALL REDUCE_DATA (NSAMPLES, COUNT, ENO ,MEAN, STO_DEV_V, STO_DEV_U, 


BLUE_GREEN) 
FINAL_U = MEAN(2 
FINAL_V = MEAN(1 
WRITE (99,¢)'Reduced velocities and std devs =’ 
WRITE (99,¢)'U dir’, FINAL_U, STO_DEV_U 
WRITE (99,¢)'V dir’, FINAL_V,STD_DEV_V 
ENOIF 
oeg@ee#ee%#etee#eeeeeee@ 
Calculate RMS values 
®m6UcaehmhCi PWUé‘C‘ Oem atwlUlcktt ete elhmUchTWmUlUcOmhUCc HUmhC ShC~™S 


OPEN(UNI T=43, FILE=’VELOCITY.DAT’ ,STATUS=’OLD’ ) 


UJ = 6.00 
VV = 6.08 
UV = 6.08 
0O 216 I = 1, END 
IF SIC recEEWAED a) READ( 43,7) eT GENL)} VORTEVEL 
IF (BLUE_GREEN.EQ.1) READ(43,7) HORIZ_VEL(1),VERT_VEL(I1) 
UU = UW + HORIZ_VEL( I )e*e2)—(FINAL_Uee2) ) 
W = W + ((VERT_VEL(I)002)-(FINAL_Vee2)) 
UV = UV + ((HORIZ_VEL(1)*VERT_VEL(1))—(FINAL_UeF INAL_V) ) 
CLOSE(43) 
UU = UW / FLOAT(END) 
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W @ W / eee 

UV = UV / FLOAT(END 

WRITE(e,e)° ' 

WRITE(¢,¢)’Reduced RMS volues’ 

WRITE(Ce,¢)'U—-RMS,V—-RMS' .UU, VV 

WRITE(e, ©)’ UV-RMS° ,UV 

WRITE(9 Me 

WRITE(99,¢)°UV=-RMS' , UV 

WRITE (25,8) X_POSITION, Z_POSITION, FINAL_U, STDO_OEV_U, FINAL_V, STD_OEV_V, 
& UU, VV ,UV 


TALLY = TALLY + 1 


Ceeeeseeeenevreeveeeeeeet ee 8 8 © © @ @ @ 
C Options to view histogrom and toke dato again 
Ceeeveev@eexsxseeeeeseeee eee & € &© 6 8 Oe 
WRITE(e,¢)' ° 
WRITE (*,¢) ‘Do you wont o Histogram? (1/8)’ 
READ (e,¢) HISTO 
IF (HISTO.EQ.1) CALL HISTOGRAM(COUNT , END, BLUE_GREEN) 
WRITE (¢,6) 
READ (¢,1) ANS 
IF (ANS.EQ.’n’.OR.ANS.EQ.°N’) GOTO 148 
GOTO 108 
148 CLOSE (25) 
CLOSE(99) 
Ceeeoeeeeeeeeeeeeeeeeeeeeee @ 
c Convert output data into proper format for 
C “PLOTCV" Vector Plot of velocity profiles 
Coeoeosoeseeeeeeeeeeeeeeee & & 6 8 Oe 
OPEN (UNTT=45, ,FILE='FINAL.OAT’ ,STATUS='OLD’ ) 
Perlite: 
OPEN( UNI T=#46,FILE='GRA.OAT’ ,STATUS=’ NEW’ ) 
WRITE(46,13) M,FREE_STREAM, ALPHA, TIME,11,16,16 
ZERO = 0.86 
DO 15@ I=1, TALLY 
READ(45,11)X_LOC(I),Z_LOC(1), UCI) .V(I),U_U(I).V_V(1).U_V(1) 
LOC UVa CZaMe Neg” 
Z_LOC( I )=(Z_LOC(1)42.5)¢.24 
vie = i / FREE_STREAM 
V(I) = V(I) / FREE STREAM 
U_U(I) = U_LU(I) / FREE_STREAM 
V v1 = Wii} / FREE_STREAM 
U_V(I) = U_V(I) / FREE_STREAM 
158 To te TURRET ae 
& U_V(I).2 
cr ece rae 
CLOSE (46 
168 WRITE (¢,5) 
G 
c Formot stotements 
c 
1 FORMAT axe 
2 FORMAT (1X,'Scanivalve § °,12) 
3 FORMAT (1X,'Port §°',12,3X, Average Voltage =°,F12.4,3X, 
& "Standord Deviation =’ ,F12.4) 
4 FORMAT (1X, Port $°',12,3X, ‘Pressure =°,F12.4,' (psi)’,3X, 
& ‘Standard Deviation =',F12.4,° (psi)’) 


5 FORMAT (//1X,'This is the ENDI’//) 

6 FORMAT (/1X,'Move taskmoster again and toke data? (y/n) [DEF - Y}'/) 

7 FORMAT (1X,2F11.4) 

8 FORMAT (1X,2F7.2,F10.2,F8.2,F10.2,F8.2, 1P3E12.4) 

9 FORMAT (11) 

11 foe Tae ee ae 

12. FORMAT(1X,13,1P10E12.4) 

13 FORMAT(1X,'GRA.DAT, MACH NO = °,F4.2,', FREE STREAM VELOCITY = ',F5.1, 
", ALPHA = ° ,F4.1,/1X, September 9, 1987’ ,4x,14,° hrs.’ 

//4%,"NPROF NSAMPL  NCOEF’/7X,12,6X,12,6X,12/1X,' ID —- X-POS", 

: Z—POS U Vv Ww wee’. 


fe fe Oe fe 


V92 wer2 UV Uw') 
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14 ee , 
16 FORMAT(6X,°x' 6X, °2',12K,°u’, 17K, *v'’, 11K, "ut2’ 9X, 'vt2", 10x, 'uv’) 

STOP 
c 
C Scanivaive Error Statements } 
c 
170 IF (HOME) THEN 

WRITE (*,*)’Problem in homing Scanivaives' 
ELSE 
WRITE (*,¢)’Problem in stepping Scanivaives’ 

ENO IF 

IF oleaist WRITE ey transfer failed in J-Board’ 

IF (1ER.€Q.3) WRITE (¢,¢)'Failure in detaching J-Board’ 

IF (HOME) GOTO 1@ 

GOTO 22 


ENO 
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f. 


SUBROUTINE SETUP(CLK_ID,DEV_ID,RATE,NCHANL NPT, NPTS) 





Subroutine written by Maj Richard R Ryles, NPS Student, 13 May 87 
Updated 15 Jul 87 


This subroutine is intended to attoch the internol clock ond AD boord 
to the main program and aliow them to be used for dota acquisition. 


This subroutine should be used with subroutines GETOATA, OETACH 





INCLUDE ‘SYS$LIBRARY:LIOSET.FOR’ 





Deciore tocol voriables | 


NOON ANAINANNDIAANAYA 





CHARACTERe4 CLK_DEV 

CHARACTERe1 ANS 

INTEGER SYS_STAT 'Stotus ret by LIO calls 
INTEGER ODEV_IO,CLK_IO 'LIO device I0 

INTEGER NCHAN ‘Channel # 

INTEGER NPT ,NPTS 

REAL RATE, TIME 





Input desired clack and somple porometers | 


QOOaND 





WRITE(¢,1) 
REAO(«,2)NCHANL 
12 WRITE(¢,3) 

READ(«,¢)RATE 

IF(RATE.GT.60@.00) THEN 
WRITE(*,4) 
GOTO 10 

ELSE 
IF(RATE.LT.0.01)WRITE(*,5) 
IF(RATE.LT.@.01)GOTO 10 

ENOIF 

20 WRITE(¢,6) 

READ(e«,7) NPT 

IF (NPT.GT.5000) THEN 
WRITE(*,8) 
GOTO 20 

ELSE 
IF(NPT.LE.O)WRITE(¢,9) 
IF(NPT.LE.@)GOTO 20 

ENOIF 

NPT = NPT + 1 

TIME = FLOAT (NPT-1)/RATE 

WRITE(*,11)TIME, TIME/68. 

NPTS = 2 ¢ NPT 

CLK_ODEV = ‘KZAQ' 

WRITE(*,12) 

READ(«,13)ANS 

IF (ANS.EQ.’y’.OR.ANS.EQ.’Y') THEN 
WRITE(*, 14) 
READ(«,15)CLK_DEV 

ENOIF 


Attach the KWV11-C Clock ond set up for mapped 1/0 | 


QAaAD 


SYS_STAT = LIO$ATTACH (CLK_IO, CLK_OEV, LIO$K_MAP) 
IF(.NOT. SYS_STAT) CALL LIBSSIGNAL(XVAL(SYS_STAT)) 
SYS_STAT = LIO$SET_R (CLK_IO, LIO$SK_CLK_RATE, 1,RATE) 
IF(.NOT. SYS_STAT) CALL LIB$SIGNAL(%VAL(SYS_STAT) ) 
SYS_STAT = LIO$SET_I (CLK_10. LIO$K_TRIG, 1, LIO$K_IMMEDIATE) 
IF(.NOT. SYS_STAT) CALL LIBSSIGNAL(%VAL(SYS_STAT)) 


OaND 


Attach the ADV11-D ond set up for mapped 1/0 | 


SYS_STAT = LIO$ATTACH (DEV_IO. ‘'AZA@’, LIO$K_MAP) 
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IF(.NOT. SYS_STAT) CALL LIBSSIGNAL(XVAL(SYS_STAT)) 


C 
C Set up the ADV for synchronous transfer | 
C Set chonne! (nchan! ) | 
C Gain of 1 | 
Cc Trigger on LIO$READ and fill buffer as fast as possibile | 
c 
SYS_STAT = LIO$SET_I (DEV_ID, LIO$K_AD_CHAN, 1,NCHANL, ) 
IF(.NOT. SYS_STAT) CALL LIB$SIGNAL(%VAL(SYS_STAT)) 
SYS_STAT = LIO$SET_I (O€V_IO, LIO$K_AD_GAIN, 1, 1,) 
IF(.NOT. SYS_STAT) CALL LIB$SIGNAL(%VAL(SYS_STAT)) 
SYS_STAT = LIO$SET_!I (DEV_ID, LIO$K_TRIG,2,LIO$K_CLK_POINT,CLK_IDO) 
IF(.NOT. SYS_STAT) CALL LIB$SSIGNAL(XVAL(SYS_STAT)) 
c 
c Format Statements 
€ 
1 FORMAT(1ix,'Enter channel number to sample: (0,1,8,9)') 
Z FORMAT (12) 
3 FORMAT(/1x,'Input desired Clock Rate (600 Hz -— 0.01 Hz):'/) 
4 FORMAT(1x, ‘Following directians! Input too high.'//20x,'Try again’/) 
5 FORMAT(1x,'Too low! Lots of time on your hands?'//2@x,'Try again'/) 
6 FORMAT(/1x,'No. of data points: (any number less than 5,001)’/) 
7 FORMAT(14) 
8 FORMAT(/1x, ‘Hey! Follow directions! Too high. Try again!'/) 
9 FORMAT(/1x,'Hey! Last input is nonsense! Try again!'/) 
1 FORMAT(/1x,'The approximate run time for each port wil! be ',f6.1, 
& ' seconds (or ',f5.2,' minutes)'/) 
pz FORMAT (/1x, ‘Do you want to change the Clock Device Name from [KZA0]?', 
& (y/n) [DEF - noJ}'/) 
13 FORMAT (A1) 
14 FORMAT(1ix, Enter Clock Device Name: (KZA@,KZB®) ') 
15 FORMAT (A4) 
C 
RETURN 
ENO 
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SUBROUTINE GETDATA(DEV_IO,NPT,NPTS, VOLTAGES) 





Subroutine written by Maj Richord R Ryles, NPS Student, 28 Apr 87 | 
Update 15 Jul 87 
| 
This subroutine is intended to acquire the dato from the A-D boord | 
with the specified porometers from subroutine SETUP. It also will 
transtate the A-O data into voltages to be converted to physica! 
quantities in the main program. | 
This subroutine should be used with subroutines SETUP,DETACH,OUTPUTO | 
INCLUDE 'SYS$LIBRARY:LIOSET. FOR’ 
Declare local variables | 
INTEGER SYS_STAT ‘Status ret by LIO calls 
INTEGER DEV_IO 'LIO device ID 
INTEGER DATA_LENGTH 'Number of bytes of dato read 
INTEGER NPT,NPTS 'Number of points to sample 
Declare data buffers | 
INTEGERe2 RAW_DATA(5@@@) 15800 point row data buffer 
REAL VOLTAGES (5820 ) '5Q@0@ voltages output 
Get o raw_doto buffer of ‘npt’ values | 
Uses LIO$READ to read the A/D | 
Note length of buffer is in bytes | 
As is the returned data length | 
SYS_STAT = LIO$READ(DEV_IO, RAW_DATA, NPTS, DATA_LENGTH, 
IF(.NOT. SYS_STAT) CALL LIB$SIGNAL(XVAL(SYS_STAT) 
Convert the raw dato to voltages using LSP$FORMAT_TRANSLATE_ADC | 
CALL LSP$FORMAT_TRANSLATE_ADC(RAW_DATA, VOLTAGES, NPT, , , ) 
RETURN 
ENO 


66 


SUBROUTINE DETACH(CLK_I0,D0EV_IO) 


Subroutine written by Maj Richord R Ryles, NPS Student 28 Apr 87 
Updote 16 jun 87 


This subroutine is intended to detach the internal clock and A-—D board 


This subroutine should be used with subroutine SETUP 


MDQAAANANANANNOA 


INCLUDE 'SYS$LIBRARY: LIOSET. FOR’ 


Ceeeesesee sess SeseSSESESSTHKEseseeseseseeeeseseeeseeeseseseseseseeseees 


G Declare local variables 
Ceeeeeee sees eS eseSeFeeSeeSeeesseseseesesseesesesesseesseresensessesesese 
INTEGER SYS_STAT 'Status ret by LIO calls 


INTEGER DEV_IO,CLK_ID 'LIO device 10 
CeeeSHee eS SSSSSSHSHSSHSSSSSSHHSHSSHSSHSKCHSHSHHSHFesesseseseaseseseuseossesse 
C Detach from the Clock 
CeSeeeS SST HSHSSSHSTHSHSHSHSHTSHHSHSTHSSHSSHHSHSHSHSHSSHKSsseseseseseeeesesessesesesece 
SYS_STAT = LIO$DETACH (CLK_ID,) 
IF(.NOT. SYS_STAT) CALL LIB$SIGNAL(%VAL(SYS_STAT)) 
Ceee sss SSSESKSHSHHEHSSHSHSSHSHSHSHSHESHSSHSHSSSeseSeesesseesesenesessesenesessess 
C Detach from the A/D 
Cees ses SSHSSSKSSTSTSSHSHSHSSHSHSHSHESHSHSHTHSHSESESEesesesesesesesesseesseseueeee 
SYS_STAT = LIO$DETACH(DEV_IDO, 
IF(.NOT. SYS_STAT) CALL LIB$SIGNAL(%VAL(SYS_STAT) ) 
WRITE(*,1) 
1 FORMAT(/1ix, ’Both clock ond A-D device are detached'/) 
RETURN 
END 
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from the main progrom and allow them to be used for others application. 





SUBROUTINE OUTPUTD(RATE.NPT .NCHANL, AMEAN ;STD_DEV, J, VOLTAGES) 





C 
C Subroutine written by Maj Richord R Ryles, NPS Student 28 Apr 87 
C Updote 24 Jul 87 
C | 
C This subroutine is intended to output diagnostic results too file. | 
C | 
C This subroutine should be used with subroutine SETUP, DETACH, GETDATA | 
g 
C Declare local variables 
G 

INTEGER NCHANL 'Channel § 

INTEGER J 'General loop indexes 

INTEGER NPT 'Number of points 

REAL «8 AMEAN ,STO_DEV 

REAL RATE, VOLTAGES( 500) 
C 
C Output dota to a file 
C 

OPEN(UNIT = 10, FILE = ‘SCANIDIAG.DAT', STATUS = 'NEW') 

WRITE(1@,1) J, RATE, NPT-1, NCHANI, AMEAN, STO_DEV 
1 FORMAT (/1ix, 'Scanivaive $',i3/1x, Clock Rote =',f6.1,5x, Dato Points’, 

& ‘Sompled =',i4/1x, Channel’ ,5x, Mean’ ,4x, Standord Deviation’/ 
& 3x,i13,2f12.4/) 

WRITE(1@,¢)'Scanivalve #',J 

DO 20 [=1,NPT 
20 WRITE(10,¢)’I=’,1,VOLTAGES(1) 

CLOSE( 10) 

RETURN 

END 
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Cessesseseseeseeesesesses 


PURPOSE : 
Orives the acqui 


ERROR HANDLING: 
IER = 1 means no "f 


EXTERNAL REFERENCES: 
DVWA Calis ORV 


DEVELOPMENT HISTORY: 
DATE INITIALS 
@7/26/83  TML 
12/12/85 CLH 


@2/25/87  GBG 
@5/06/87  CLH 
@6/08/87 RRR 
@7/15/87 RRR 


@7/24/87 RRR 
08/29/87 RRR 


DQOAIAIAIIAIAIAIANAAANAANADIANAIANAANANANAANDNAN Oo 


C AUTHOR(S): Ted Lichten 


C—eeeseeeeeeesseseececes 


C 
IMPLICIT NONE 


SUBROUTINE HWWA(NSAMPLES , ANSWER , BLUE_GREEN) 


sition of the LOV data. 


atai" error. 


NAME DESCRIPTION AND SOURCE 


11-WA device driver to acquire A/D data. 


DESCRIPTION 
Adapted from DAPNT 
Modified to print results in either 
octal or integer. 
Modified for use with ORV11-WA and VWS 
Modified to accept "new" DAMUX1. 
Modified to subroutine + deleted CLH 5/87 Mod 
Modified to calculate velocities + optimized 
for speed of execution 
Modified to view Histogram 
Modified to accept calibration constant 


stein, Informatics General Corp. 
Seesoeeoeeeseseeeeeeoseeeeeseeeeseeese eee eeeeeeeea ee ee eee se 


c 
PARAMETER NRAWMX=8 0000 
PARAMETER LUNT [25 
PARAMETER LUNOUT=6 
g 
INTEGER*2 IRAWDT (NRAWMX ) 
INTEGER®2 11RAWDT (NRAWMX) , I2RAWDT (NRAWMX ) 
INTEGER IER, HISTO,Nf r,ONE, TWO 
INTEGER ITYESNO, NRAWIN, I, 8BLUE_GREEN 
INTEGER NCHANS, NSAMPLES, INDEX, JNDEX 
INTEGER NS . ANSWER , COUNT , END 
REAL TIME (NRAWMX) , FREQ(NRAWMX) , VELOC (NRAWMX ) , I 3RAWDT (NRAWMX ) 
REAL MIX(2),GAM(2) ,DEL(2) ,KAP(2) ,0(2),1(2) ,MEAN(2) , FREQUENCY(2 , 8001) 
REAL FREQONE , FREQTWO, FRE(NRAWMX) , FREQ1_AVE, FREQ2_AVE, CAL_CONST(2) 
LOGICAL DATSELCT (3) 
CHARACTER#19 DIRECTION 
CHARACTERe 1 ANS 
c 
DATA MIX(1).MIX(2) / —1.E-6,-1.E-6/ 
c 
INCLUDE *($SYSSRVNAM) ’ 
c 
1800 CONTINUE 
Ceeeeeeeeteeeeeeenstse ee ee ee ehh hthlchlUlcrhlUl lt 
C Fill the raw data buffer with this bit pattern: 
C *180101010180101010" 
Cc This equals 125252(8), -21846(10). 
Ceee ewe eeeeeeeeeeeeeeeeree es eee # 86 
DO 200 | = 1, NRAWMQX 
IRAWOT(I) = '125252'0 
11RAWOT(I) = '@@00@0'0 
I2RAWDOT(I) = °'@0@@00'o 
13RAWDT(I) = @.@ 


20@ CONTINUE 


QOaM 


Input number of samples to be acquired 
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NCHANS = 4 
300 WRITE (LUNTI, 18) 
12 FORMAT (' Enter number of samples:’) 
READ (LUNTI. «) NSAMPLES 
NRAWIN = NSAMPLES « NCHANS 
IF (NRAWIN .GT. 32767 .OR. NRAWIN .LE. @) THEN 
WRITE (LUNTI, 20) 


220 FORMAT (’ The number of samples « number of channels must be 32,767’ 
+ "or less.') 
GO TO 300 
ENO IF 
Ceeeeeet ee eeseeeeeweeeeeeeeeee 8 @ @ 
C Set up the LOV parameters prior to acquisition 


G | 
IF (BLUE_GREEN. EQ. 1) THEN 
GAM(1) = @.488 
GAM(2) = @.5145 
IF(CAL_CONST(1).EQ.1.) THEN 
WRITE(e,*)' Enter Blue calibration constant’ 
READ(«,¢) CAL_CONST(1) 
ENOIF 
IF(CAL_CONST(2).£Q.1.) THEN 
WRITE(e,¢)'Enter Green calibration constant’ 
READ(*,*) CAL_CONST(2) 
ENDIF 
ELSE 
GAM(1) = @.5145 
GAM(2) = @.488 
IF(CAL_CONST(1).EQ.1.) THEN 
WRITE(*,¢)'Enter Green calibration constant’ 
READ(*,*) CAL_CONST(1) 
ENDIF 
IF (CAL_CONST(2).£Q.1.) THEN 
WRITE(e,*)'Enter Blue calibration constant’ 
READ(*,*) CAL_CONST(2) 
ENOIF 
ENDIF 
IF(MIX(1).LT.0.0) THEN 
IF(BLUE_GREEN.EQ.1)WRITE(*,¢)’ Input BLUE mixing frequency (MHz) ' 
IF (BLUE_GREEN.EQ.@)WRITE(*,*)’ Input GREEN mixing frequency (MHz) ' 
READ(*.*) MIX(1) 
ENOIF 
IF(MIX(2) LT.@.0) THEN 
IF(BLUE_GREEN.£Q.1)WRITE(«,*)’ Input GREEN mixing frequency (MHz)’ 
IF(BLUE_GREEN. £Q.@)WRITE(e¢,¢)'Input BLUE mixing frequency (MHz)' 
READ(e,*) MIX(2) 
ENOTF 
IF(ANSWER £Q.@) THEN 
WRITE(*,*)'Do you want probe volume specifications? (y/n)' 
REAO(e.’(a)*) ANS 
IF(ANS EQ. ‘'Y'.OR.ANS.£Q.‘y’) ANSWER = 1 


ENO IF 
WRITE(*,¢)' 
C 
C Select the data to be returned. 
c 
DATSELCT (1) = FALSE. 
OATSELCT (2) = FALSE. 
OATSELCT (3) = .TRUE. 
Ces et eee ee eee eeeeee ee eel ehl hl OhUlUc MhlUc hl lt 
C LOV Dato Acquisition 
Ceee ee © ee we ee ee hl hl OhUlUlcrOhUc OChUlUc OChUlcr OmhUc OhUlUc OhlUcr OhlUhrOMhlUh Hh Uhh 
ONE = 1 
TWO = 2 


OPEN(UNIT#=5,FILE=’VELOCITY.OAT’ ,STATUS='NEW' ) 

ODO INDEX = 1, NSAMPLES 
CALL DAMUX (NCHANS, ONE, DATSELCT, , . IRAWDOT, IER) 
IF(IER .NE. 1) WRITE(*,¢)’ TER in OAMUX =’, TER 


TO 





Cae ele SueMNeIGIGIStS GS 5 a6 6 © 0 6 6 6 0 8 8 ols 8 


C Conversion to Fraquency 
Cesteoeeerteoeoeeeeeeeeeeeeteeeeseeeeeee 
0O [=1,2 


IRAWOT(I) = INOT(IRAWOT(I)) 
IZSRAWOT(I) = 0.2 


J2RAWOT(I) = @ 
]1RAWOT(1) = @ 
]1RAWOT(1) = IIBITS(IRAWOT(I),10,4) 
I2RAWOT(1) = LIBITS(IRAWOT(I),@, 10) 


I3RAWOT(1) = (2.eeI1RAWDT(I)) ©( float(I12RAWDT(I))) 
TIME(1) = I3RAWOT(1)/32. 
FREQ(I) = 1000./TIME(1) 
END DO 
FREQUENCY(1, INDEX) = ae 
FREQUENCY(2, INDEX) = FREQ(2 
es e@ee#ee#ee#e##ese#e8ee08e#e88 #808 080880 80800hCUmM8t6UmUC<“SUCUMDSLUC UCU SUCUMOUUCUD 
Canvert frequency data to velacities 
oeost 8@ oe s8kts 8 8s: &®.lmUcCOrhC<CNSHDMC N H]WDMCU(CW ]hMhM—COCODUCUCc TUc Uc OrhUC TU OO TChUhOOmUCUc HOTmUC OrhUChCrFmhlhUD 
CALL VELOCITY (FREQ, VELOC, TWO,MIX, ANSWER ,GAM,DEL,KAP,D,1,Nfr,CAL_CONST ) 
WRITE(S,30) VELOC(1),VELOC(2) 
38 FORMAT (1X,2F11.4) 
END DO 
CLOSE(5) 
ees e0es8s:e8 e0 8¢ 8 800 8 00 8080808080880 80ChlUmM0—606—C<=UUCUOUmUmCOUCUC UC OUC UC OUCUCUCOUCUmCOUmUC HDLC CUD 
Calculate average caunter frequency and print to sceen for diagnostics 
a a 
FREQONE = 0.0 
FREQTWO = 0.0 
DO 34 | = 1,NSAMPLES 
FREQONE = FREQONE + FREQUENCY(1,1) 
34 FREQTWO = FREQTWO + FREQUENCY(2,1) 
FREQ1_AVE = FREQONE/FLOAT(NSAMPLES ) 
FREQ2_AVE = FREQTWO/FLOAT(NSAMPLES) 
WRITE(e,¢)’ ' 
WRITE(*,*) Average counter frequency § 1 
WRITE(*,*)’Average counter frequency § 2 


OaND 
e 


OOO 
. 
. 


=', FREQI_AVE 
=', FREQ2_AVE 


Ce eeseerteees#e#e eee ee © @ e@ © &@hchmhUcrMrmhCUCc OhUCUCOhmhUCUCcrhUCUchUlUcOhUlUcrhUlUcrhUlUh hl 
G Write Probe Volume dimensions and porometers toa file 
Ceeererteveeeeeteeeeeeseeee sts © © © & & & oo 
IF (ANSWER. EQ. 1) THEN 
OPEN(UNI T=15,FILE=’Probe. dat’ ,STATUS='NEW’ ) 
WRITE(15,°)'Probe Volume Data’ 
WRITE( 1S, ©) 0 —_— 
& ee 
DO 40 I = 1,2 
IF(BLUE_GREEN.EQ.1) DIRECTION = ‘Horizontal! Velocity’ 
IF(BLUE_GREEN.EQ.@) DIRECTION = ‘Vertical Velocity ' 
IF(I.€Q.2) THEN 
IF(BLUE_GREEN.EQ.1) DIRECTION = ‘Vertical Velocity 
IF (BLUE_GREEN.EQ.@) DIRECTION = ‘Horizontal Velocity’ 
ENDIF 
WRITE(15,1)GAM(1) ,OIRECTION,DEL(1),KAP(I),O(1).!(1)*.001.Nfr 
1 FORMAT(/1X, 'Wavelength af light <microns> = ',f7.4,' (Measuring’, 
& a19,')'//1x,'Fringe Spacing <microns> = ',f10.2//ix,'Lens’, 
& angle <deg> =',f7.2//1x,'Probe volume diameter <microns> = ', 
& f10.5//ix,’ Length of measuring volume <mm> = ',f10.4//1x, 
& ‘Number of Fringes =',i3/) 
WRITE( 15,8)’ 
& ® 6 
40 CONT INUE 
CLOSE( 15) 
ENDIF 
Cs s ss ss eos ee se emhlUcrhmhlUc OhUMMChUcUOOUC HOCUhMH HP FP eS Fe Fe oe ew se Fh Fe Fe Fe Fe Fe ee Fe oO Fs 
C Option ta print a histogram of velocity data an screen for diagnostics 
Ceeseseseee#es8see#e#ee#eeses#seeeees © © © © & & & Oe Fe Fe 8 
WRITE(*,#)’ ° 


WRITE(*,*)’Do you want a Histagram? (1/0)’ 


A 


A000 


READ(«,*) HISTO 
END = NSAMPLES 


IF (HISTO.EQ.1)CALL HISTOGRAM(COUNT , END .BLUE_GREEN) 


Again? 


WRITE(LUNTI,*)’ 

WRITE (LUNTI,*) ‘Acquire LDV Data again? 
READ (LUNTI,*) IYESNO 

WRITE(*,¢)’ ' 

IF ( IYESNO.EQ.1 ) GO TO 102 


RETURN 
END 


(O=NO , 


2 


imYES) 





NQANANAAANAND 


SUBROUTINE VELOCITY( FREQ, VELOC .N,MIX, ANSWER ,GAM,OEL,KAP.De. |! Nfr. 
& 


CAL_CONST) 


Subroutine written by Moj Richard R Ryles, NPS Student, 9 Jun 87 
Last Update 29 Aug 87 


This subroutine was written to convert the LDV raw data from frequency 
to velocities. 


REAL  THETA,GAMMA,DELTA,d,f,E.KAPPA,Dsubm, !m,De,FREQ(N) ,VELOC(N) ,GAM(2) 
REAL SHIFT,MIX(2),DEL(2),CAL_CONST(2),KAP(2),Om(2),1(2) 
INTEGER IANS, ANSWER ,Nfr 
be erreeReets§ coueme> CAL_CONST(1) = 1.0 
IF (CAL_CONST(2).€Q.@.0@) CAL_CONST(2) = 1.0 
PI] = ASIN(1.0)e2. 
d= 13. |! Beam spacing (mm) 
f = 600. ! Focal length of transmitting lens (mn) 
| 
! 


E = 1. Beam expansion ratio 
KAPPA = ATAND(E*ed/2./f) Lens angle (deg) 
DO 5 1 #=# 1,2 


GAMMA = GAM(1) ! Wavelength of Laser (microns) 
DELTA = CAL_CONST(I) © GAMMA / (2.¢SINDO(KAPPA)) ! Fringe spacing 


De = 1.1 ! Diameter of the incident laser beam (mn) 
Osubm = (4. eGAMMAs{/(PleDeeE)) ! Probe volume diameter (microns) 
Im = Dsubm / TANO(KAPPA) ! Probe volume length (microns) 


Nfr = 1.27 © d / De 
DEL(I) = DELTA 


! Number of fringes 
KAP(1) = KAPPA ! 

! 

! 


Variables passed to output file 
Om(1) = Osubm - : a : : 
1(1) = In 
SHIFT = 48.2 
FREQ(I) = FREQ(I) — SHIFT + MIX(1) 
VELOC(i) = DELTA ¢« FREQ(i) 
RETURN 
ENO 


oe oe + os 1 


= Bragg shift (MHz) 

MIX = Mixing Freq (MHz) 

FREQ = Adjusted Doppler Freq 
(MHz ) 
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SUBROUTINE HISTOGRAM( COUNTER, END, BLUE_GREEN) 








C 
C Subroutine written by Maj Richard R Ryles. NPS Student, 16 Jun 87 
; Update 29 Aug 87 
C This subroutine is intended to create two histograms of velocity 
c data from a two component LDV. 
Cc 
Cc 
CHARACTERe4 STRING 
CHARACTERe22 LABELX ! variables for the axis 
CHARACTER?17 LABELY | tables 
CHARACTER#29 LABE LM 
INTEGER STATUS, I, ISHADE , COUNTER, END 
INTEGER N(3),BLUE_GREEN 
REAL CONDE sea a 
REAL YCONTROL( 4) , COUNT (8) 
REAL XLOw(8) 
REAL XHIGH(8) 
REAL Y(8) , VELOCITY(82@@, 2) ,x(8200) 
g 
N(1) = 1 [| number of data points 
N(2) # 1 | number of columns of data 
N(3) = 1 | number of rows of data 
Ce eoeetee54useet*@eeseseseeeeee 
C Read Velocity Data File 
C us e*eeoe<4useoeteeee ee ¢ @ @ 6 
OPEN( UNI T=5,FILE=’VELOCITY.DAT’ ,STATUS='OLD’ ) 
IF (BLUE_GREEN. EQ. 1) THEN 
DO 10 I = 1,ENOD 
18 READ(5,1) VELOCITY(1,1),VELOCITY(1,2) 
ELSE 
DO 2@ I=#1,ENO 
20 READ(5,1) VELOCITY(1,2),VELOCITY(I,1) 
ENO IF 
Ce e*e2eees8e eee eee &¢ & @ 
c Find the Min and Mox values 
Ce e*0@088e8e08 8? @¢ @ 8@hUcHhUcChCUhHhUhHMhUhShUO 
CLOSE(5) 
00 86 J #= 1,2 
DO 38 K = 1,end 
30 X(K) = VELOCITY (K,J) 
XMIN = 1000.6 
XMAX = ~-1006.0 
DO 40 I = 1,end 
1a oe a XMIN = nt 
40 IF (XMAX. LT. XC] XMAX = X(/ 
Cre ee ee ee ee ee er er er er ee er ee ee er ee | 
C Initialize X-Axis controls 
Ce eeeeee eee ee © © &€ @ @ 
XCONTROL(1) = 5.5 ! tength of x-axis, in inches 
XCONTROL(2) = XMIN ' teft most value on x-oxis 
XCONTROL(3) = XMAX ! right most value on x-oxis 
XCONTROL(4) = (XMAX-XMIN)/8. | increment of plotted points 
Ce e*e0see#e8e8e8ee#e#e#8e08 8 86 @ 
C Find y data from somple array 
Ci *“eeenskeestetee@ 8 8 &@ 6 @ 
0O 5$®@ 12 1,8 
Y(I) = 6.0 
XLOW( 1) = XMIN + FLOAT(I-1)*XCONTROL(4) 
IF(1.€Q.8) THEN 
XHIGH(8) = XMAX 
XHIGH(I)= XMIN + FLOAT(1)*XCONTROL(4) 
5e COUNT(I) = XMIN + FLOAT(I1)¢XCONTROL(4) 
C 


XTOTAL = 0.6 
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DO 6@ I = 1,end 
XTOTAL = X(I) + XTOTAL 


IF(X(1).GE. MIN. ne X(1).LT.COUNT(1)) Y(1) = Y(1) + 1. 
IF(X(1).GE. COUNT .ANO .X(1I).LT.COUNT(2)) Y(2) = Y(2) + 1.6 
IF(X(1).GE. COUNT 2 .ANO.X(I -LT.COUNT 3)) Y(3) = Y(3) + 1.0 
IF(X(1).GE.COUNT(3).ANO.X(I).LT.COUNT(4)) Y(4) = Y(4) + 1.0 
IF(X(1I).GE.COUNT (4) .ANO.X(I).LT.COUNT(S)) Y(S) = Y(5) + 1.0 
IF(X(1).GE.COUNT(5).ANO.X(1).LT.COUNT(6)) Y(6) = Y(6) + 1.8 
IF(X(1).GE.COUNT(6).ANO.X(I).LT.COUNT(7)) Y(7) = Y(7) + 1.8 
ee tuee. .COUNT(7).ANO.X(1).LT.COUNT(8)) Y(8) = Y(8) + 1.0 
60 1F(X(1).GE.COUNT(8).ANO.X(1I).LT.xXMAX) Y(8) = Y(8) + 1.0 
XMEAN = XTOTAL/FLOAT( END) 
c 
YMAX = 0.0 
00 70 I = 1,8 
Y(1) = Y(1)/FLOAT (END) 
IF (YMAX.LE.Y(1))YMAX = Y(I) 
78 Y(I) = 100. © Y(I) 
WRITE(*,2)YMAXKe100. , XMEAN 
Ces seseeseveeeseeese 8 &@ & 8 
c Initialize Y axis controls 
Ceeevevseieseeeseee8 © @ @ 
Meotrcotts = 4.5 | height of graph in inches 
YCONTROL(2) = 0.8 | lower boundary value 
i a = YMAXe120. | upper boundary value 
YCONTROL(4) = 5.0 | increment of plotted points 


Ceeewesseeesesesees es 8 ee 8 &@ @ 8 
C Histogrom labels 
Ceeevvses tees eeese ees es @ 8 
LABELX = ° VELOCITIES (mps) : ! x-axis label 
IF (COUNTER.GT.@)LABELX = ‘FINAL VELOCITIES (mps)’ | x-axis label 
LABELY = ‘PERCENTAGE OF PTS’ ! y-axis label 
LABELM = °HORIZONTAL VELOCITY HISTOGRAM’ ! upper label 
IF(J.EQ.2) LABELM = ° VERTICAL VELOCITY HISTOGRAM’ 
STRING = °IXSY’ ! mode string 
ISHADE = 2 | filling of graph 


oeeseeee#28ee2em608kmtmUmUm OU UODUCUCUHUCUC HUC TWUCUh]HOTChUCUC OrhUC OUCUhTCUC(UhHhU(Uh!HhMUC(<~h SHhU!B 
Set up coordinate system with values 
eo @eeee#e#8ee08e88ee0@8@#e0e08e88f8080806h6mM06C~<CGCi<C HUCUCLUCUCUD 
IF (XMIN.GE.XMAX)GOTO 90 
CALL LGP$PLOT(1,STRING, .0.N, LABELX, LABELY,STATUS,., 
= XCONTROL , YCONTROL, , LABELM) 
IF (STATUS .NE. LGP$_SUCCESS) THEN 
CALL LIB$SIGNAL(XVAL(STATUS) ) 
ENO IF 
eeeee#2ee#ee#e#e#e#ee#eee#e#eeseeseeeeseesstes 
Call routine and plot histogram 
@eee#2eee#e#s##e#eestes#se#ee#see##e#e @ 
CALL LGPS$HIST(1,XLOW, XHIGH, Y,8,STATUS, , ISHADE) 
IF (STATUS .NE. LGP$_SUCCESS) THEN 
CALL LIB$SIGNAL(XVAL(STATUS) ) 


Oa 


Onan 


ENO IF 
PAUSE ‘Waiting for viewing/printing Histogram; Type CONT when finished’ 
CALL LGP$TERMINATE_PLOT(1 | terminate plotting 
80 CONT I NUE 
RETURN 
99 i 2 
WRITE(e¢,¢)’eee No Histogrom created: xmindexhigh eee’ 
GOTO 82 
Ceeeveeeses eevsxst @ @ &@ © @ &@ @ @ 
Cc Formats 
Ceeeevevseeeeseeet &@ @ 8 
1 SN Liban 
2 FORMAT(/1X, ‘MAXIMUM PERCENTAGE OF SECTORS = °,F6.1//1X, ‘SAMPLE MEAN = °, 


& F9.2) 
ENO 


79 


SUBROUTINE REDUCE_DATA(NSAMP , COUNTER , ENO ,MEAN, STD_DEV1,STD_DEV2. 
& BLUE_GREEN) 





Subrautine written by Maj Richard R Ryles, NPS Student, 12 Aug 87 
Update 29 Aug 87 


This subrautine is excludes all velacity data from the LDV that fall 
outside twa standord deviatians from the mean. 


OOOaAaN00D 





REAL*®8 STO_DEV1,STD_DEV2,NEW_STD(2) , TOTAL(2),SUM(2),X.Y 
REAL VELOC(2,1@@@@) ,MEAN(2) ,UP,DOWN, TOLERANCE 
REAL NEW(2),VELO(2 , 20000) 
INTEGER NSAMP ,COUNT(2), TOT, END,BLUE_GREEN 
Ce eeeeke#weee Cteeeeeee ewe ehlcehmlUchmhUch MmhUCUcOMmhUCUCcOMrmhUCcOrhCUcrhUcUOOmhUchOUh HOU hUhCOlhlUO 
C Read velacity data from file 
Ces ese eeeeeteeeeeeeseeeesteeseeeeest © © & & 
OPEN (UNI T=3, FILE=’VELOCITY.DAT’ ,STATUS='OLD’ ) 
DO 1@ I = 1,NSAMP 
1 READ(3,1)VELOC(1,1),VELOC(2,1) 
CLOSE(3) 
es «eee8ee¢e¢6e e860 6¢e 8@ 0ee0ee80e 8 @ 8@ 60 60 6@ 60 @ 0 60 6 0@6¢ 860 060 860 860 86 86 
Calculate range af acceptable values i.e. all velocity 
data that falls within two standard deviatians af the mean 
eS SS Se Se ee | 
DO 30 J = 1,2 
IF(J.€0.1)THEN 
TOLERANCE = STO_DEV1e2. 
UP = MEAN(1) + TOLERANCE 
DOWN = MEAN(1) — TOLERANCE 
ENDIF 
IF(J.€Q.2) THEN 
TOLERANCE = STD_DEV2e2. 
UP = MEAN(2) + TOLERANCE 
DOWN = MEAN(2) — TOLERANCE 
ENDIF 
* e@eees#seF#ee¢e8e8e @ ee ee e# 0e 60 80 8e@ 8 860 @ 0 0 0 860 86 86 86 6 @ 
Flag each paint autside of twa standard deviations 
eo «ee 8 @ 8 @ @ ®@ @ @ @ @ 68@mUmUOMmUmCOMUCSUCUUC HTC OMhC DMC RHhMLUC OhUhUCFMhMhC<CSTMUCTHLhUCUC OhUCChlUr 
DO 2@ | = 1,NSAMP 
IF(VELOC(J,1).LT.DOWN.OR.VELOC(J,1I).GT.UP)VELOC(J.1I)=1000.0 


eoekek® 


eg & 


CONTINUE 


Calculate new mean and standard deviatians 
Cr ee | 


OAaAAWHNH aan 


® e e e e es 
DO 60 J = 1, 
TOTAL(J) 
COUNT(J) = @ 
SUM(J) = @.De 
DO 50 1 = 1,NSAMP 
IF(VELOC(1,1).NE.100@.0.AND.VELOC(2,1).NE.1000.@) THEN 
Pons = TOTAL(J) + VELOC(J,1) 
COUNT(J) = COUNT(J) + 1 
ENDIF 
52 CONT INUE 
IF (COUNT(J).£Q.@) COUNT(J) = 1 
NEW(J) = TOTAL(J)/FLOAT(COUNT (J) ) 
DO 55 K = 1,NSAMP 
IF(VELOC(1,K).NE.10@0.0.AND.VELOC(2,K).NE.100@.@) THEN 
X = VELOC(U,K) — NEW(J) 
Y = XeXe(1./FLOAT(NSAMP-1) ) 
SUM(J) = SUM(J) + Y 
END IF 
55 CONT INUE 
NEW_STO(J) = SORT(SUM(J)) 
Q CONT INUE 


6 
c 
C Output new values to screen 
C 


®.De 


Oh a hr 


eeke 


O00 


828 


180 


N—- O00 


a eee 
WRITE(*,*) ‘New cleaned—up data’ 
IF (BLUE_GREEN.E€Q.1) THEN 
WRITE(*,¢)°U dir’ ,NEW(1),NEW_STD(1) 
WRITE(*,*)°V dir’ ,NEW(2),NEW_STD(2) 
ELSE 
WRITE(*,*)'U dir’ ,NEW(2) .NEW_STD(2) 
WRITE(*,*)°V dir’ ,NEW(1),.NEW_STD(1) 


ENDIF 
MEAN(1) = NEW(1) 
MEAN(2) = NEW(2) 


STD_DEV1 = NEW_STD(1) 
STO_DEV2 = NEW_STD(2) 
* ¢© © eh elhCUctmhC HhUC HhUC<C hUC HhmhUCr HhUC HhUC HH hUCU HMhUC HhUC HDhUC HMhC<C TCU HhUCch HhUCCmHhCOrmhHrhUCc H]hUC HOMhUC HhUC<C SMC h—C<C F!hlUh 
Write interim data file containing just good data 
Sr ee er ee Se 
OPEN(UNIT=5,FILE=’INTERIM.DAT’ ,STATUS=' NEW’ ) 
DO 70 J = 1,2 
DO 7@ I = 1, NSAMP 
IF(VELOC(1,1).NE.100@.0.AND.VELOC(2,1).NE.10@@.@) THEN 
WRITE(S,2)VELOC(J,1) 
ENDIF 
CONTINUE 
CLOSE(S) 
Cr ee 
Write final velocity data file for plotting 
* © @e8#e8¢e¢#e0e6mUMtmC UC UC HEUC HUC TD ThUCh HhhUCUc HTmhUC HMC HmhC<C HO HTUCU HhUC HOmhUCUMHMhmhCc HOhmhUC HMhUC MhC—C HMC HhMhCh—é~C<“‘ SSCP 
OPEN(UNIT=7,FILE=’VELOCITY.DAT’ ,STATUS=’NEW’ ) 
OPEN(UNIT=9 ,FILE='INTERIM.DAT’ ,STATUS#’OLD’ ) 
TOT = COUNT(1) + COUNT(2) 
DO 88 | = 1, TOT 
IF(I.LE.COUNT(1) ) THEN 
READ(9,2) VELO(1,1) 
ELSE 
READ(9,2) VELO(2,1) 
ENDIF 
CONT INUE 
CLOSE(9) 
END = COUNT(1) 
COUNTER = 1 
DO 10@ [| = 1, END 
WRITE(7,1) VELO(1,1),VELO(2,COUNT(1)+I) 
CLOSE(7) 
RETURN 


Format statements 

Cr rr ee a ee | 
FORMAT (1X,2F11.4) 

FORMAT (1X,F 11.4) 

END 


rg 
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